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ABSTRACT 


The radiational model employed here computes the planetary 
albedo and the solar absorption by atmospheric layers and by 
earth's surface for the primitive equation model of FNWC. 
meree-sScale cloud parameterization in several layers enters 
into these computations. The solar computations are made on 
a gridpoint basis employing the water-vapor mass over each 
gridpoint. Longwave cooling effects over the same layers and 
at the earth's surface are calculated using emissivity formulas 
after Sasamori, and are also dependent upon the cloud- 
parameterizations. 

Two forms of the cloud-parameterizations were tested using 
FomG gridpoint @ceanic data for 16 January 1974. The objective 
was to determine the parameterization which better verified 
the radiational balance as a function of latitude at the 
tropopause, as compared with the January satellite climatology. 
Separate heat-budget checks of the tropospheric columns and of 
the oceanic-water mass were tested for consistency with January 
Climatology. The smaller cloud-parameterization values gave 
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I. INTRODUCTION 


This thesis represents a continuing study of a radiative 
heating parameterization for use in appropriate o-layers of 
the Fleet Numerical Weather Central (FNWC) numerical predic- 
tion system. The immediate goal of this thesis is to analyze 
the heat balance of the earth-ocean system utilizing FNWC 
eumaded data-fields at constant pressure levels for 16 
January 1974. These gridpoint data were selected along four 
meridians, three in the Pacific and one in the Atlantic 
Weean. Most ot the gridpoints were located in the Northern 
Hemisphere (Fig. 1). 

Mie most important influence on the radiational disposi- 
tions (short and long-wave) of this study was the specifica- 
tion of amounts of clouds in two specific layers. The 
initial specification of the fractional amounts CL(1) and 
CL(2) are based on large-scale formulations developed by 
Smagorinsky [1960] and first used in a similar study based 
upon 16 October 1973 data by Warner [1974]. 

The governing equations of the radiational transfers of 
the model have been devised by Martin [1972, 1974], who 
modified the streams of radiation, both solar and terrestrial, 
to respond to the presence of clouds in specified amounts in 
the two layers. The radiational model developed by Martin 
milcewtilized here has similarities to those in use in UCLA 
and NCAR General Circulation models, but pane ce el adapta- 


tion to the o-levels of interest to FNWC. 


ay, 





(1,463) 


(9,55) 






Aw 145°E 


3 
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Peeure 1. FNWC polar stereographic grid and meridians (lines 
mm 2, 3, and 4) selected for study. The longitudes A are 


Shown for each meridian as well as the extent considered of 
each meridian. 
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Wormer ound thateies global albedo estimates were too 
high as compared with the satellite climatology of Raschke 
Eee act. [1973], thus leading to underestimates of the radia- 
mive balance both at the top of the atmosphere and at the 
surface for the data of his experiment. In this study, the 
original Smagorinsky [1960] cloud-specifications CL; and CL, 
were computed for use in the model; however, reduced cloud- 
amounts (2/3-CL) in each layer were also included for radia- 
tive budget comparisons. 

The comparative results afforded by the "2/3-CL" model 
gave reasonably close agreement as a function of latitude 
Wath the radiative climatology of Raschke et atl. [1973], 
fomecne NIMBUS III period 21 January - 3 February 1970, 
Svamined by these experimenters. 

In order to analyze the effects of the '"2/3-CL" versus 
fem tul)-CL model used by Warner, oceanic and tropospheric- 
Wear Dalances were also computed for both sets of cloud con- 
ditions for 16 January 1974. It was necessary in this 
Sommection to augment the purely radiative model by a turbu- 
lent boundary-layer model for sensible and latent heat 
transports as adapted from Kaitala [1974]. 

In the 2/3-CL case, the zonally averaged results for 
the components of the ocean-tropospheric system gave realis- 
tic results when compared to climatology for this mid-winter 
meriod. 

The possible error-source of using one set of computation- 
meecdata to be representative of mid-January climatology is 


recognized. Nevertheless, the heat-package results computed 


el 





with the 2/3-CL cloud model were encouraging when examined 
meainst climatology. The major effect noted by the reduction 
in cloud cover was the reduction in global albedo in both 
low- and mid-latitudes thus contributing to greater surface 
net heating rates. This latter result is in general accord 
with recent observations of Von der Haar and Hanson [1969], 
Von der Haar and Oort [1973], as well as with Raschke et al. 
[1973]. However, the specification of layered amounts of 
Clouds in the tropics still seems to give somewhat excessive 
global albedos. It is felt that the geometry of solar radia- 
tion streams impinging upon cumulus subgrid cells in the 
tropics must be reexamined so as to divert a greater per- 
centage of solar radiation downward, as compared with cor- 
responding cloud-layering reflective effects in middle and 


high latitudes. 
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II. DATA PREPARATION 


Temperature and humidity data utilized in this study 
were taken from gridpoint soundings along four meridians 
of the FNWC 63-by-63 gridmesh. The data gridpoints along 
these meridians (Fig. 1) were chosen so as to ensure condi- 
eons representative of oceanic locations for computations 
of the radiational and other heat budget items. The reasons 
for choosing oceanic locations for such computations were 
threefold: (1) the maritime-area heating-rate computations 
are likely to be relatively stable timewise over the data- 
faye lGo January 1974): (2) the constant o = P/n surfaces 
of the FNWC primitive equation system are close to being 
Semlstant pressure levels; and (3) values of the heat-budget 
items computed here may be compared with corresponding cli- 
matology during mid-winter. Corresponding values over land 
areas appear to be less amenable to comparison. 

mone the £our oceanic meridians used, a total of 67 
gridpoint soundings were selected from I = 1 to I = 25 on 
imes 1, 2, and 3 (over the Pacific Ocean), and 26 soundings 
from I = 63 to I = 38 along the Atlantic meridian, line 4 
(ere, Fig. 1). The device of selecting gridpoint soundings 
chosen from the indicated meridians of the FNWC polar 
Stereographic map, made it unnecessary to resort to spatial 
agiterpolation between original data gridpoints along the 
meridians. 

The gridpoint soundings from the data lines just mentioned 


were taken from the original FNWC 63-by-63 gridpoint analyses 


25 





Of temperature at the surface and of T(P) and q(p) at nine 
other standard pressure levels up to and including 100 mil- 
libars (Table I (a)). Those soundings were then modeled 
into a radiative sounding model compatible in the vertical 
men ENWC primitive equation prediction levels (Fig. 2). 

The Pacific Ocean gridpoint data cases were taken from 
the OOOOGMT, 16 January, 1974 FNWC analyses, while the 
Atlantic Ocean data cases were taken from the 1200GMT, 16 
wemvary $1974 analyses. The reason for the choice of the 
QOOOGMT time-analysis was to correspond to a set of "sound- 
ings'' as close as possible to local solar noon over the mid- 
Pacific, whereas those of 1200GMT approximate a near-noon 
conditions over the Atlantic meridian. 

mAeweach eridpoine selected, the original data listed the 
five moisture values given dewpoint depressions for the 
aatysis Levels from 925 mb to 400 mb. Since the standard 


micenument level vapor pressure (e ) was missing from the 


assis 
original FNWC analysis at all gridpoints, Caiy was app roo 
mated using the FNWC field of e,, a computed value of vapor 
pressure at about 20 meters above MSL in the turbulent 
boundary layer. The FNWC values of ie and e, over the ocean 
had been made available using an operational planetary bound- 
ary layer model detailed by Kaitala [1974]. 

In Sections V and VI it was necessary to recompute values 
Oe ey (and of i) by application of the FNWC planetary bound- 
ary layer model to data reported in the gridpoint soundings. 


It was found that FNWC values of e, as Hitter ktyerenonted were 


PeGurate Enough to be representative of the near-surface 
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identifers U(0), C0) O-level 
Se GO, 
Layer (0.2) a Voie ae 01 
7 nr Ge2E) en 2) G02 
(2,4) aa Se mq == = 0.3 7 
| EC) ey 
(4.6 eae | 
- @8) y--------- qr --- 7-07 
8 AGE eee) 0-08 
CL (2) 
(8 i0) PE ae 09 
| oy ene) mC O NE SES 1G 
Figure 2. Five-layer radiative sounding used in this study. 


Levels are identified by their values on the k-scale, while 
layers are identified by their level boundary indices in 
parentheses, e.g. (8,10). Pressure-Scaled water vapor and 
CO. mass increments U and C, respectively, are integrated 
with respect to the surface and introduced at even levels 
while the mixing ratio, q, is formulated at odd levels. The 
temperature, T, is expressed for all levels. Amounts of 
Clouds CL(1) and CL(2) in the layers shown have been param- 
Peer zed fOr Consideration of their radiative effects. 
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Table I(a). Example of typical FNWC sounding for gridpoint 
(1,1). The humidity parameters between 925,...,400mb are 
dew-point: depressions, and at the surface and top of the 


constant flux layer (9999*) are vapor pressures. 


Pressure (mb) TGC ) Humidity Parameters 
surface, 1000mb PAS ee: e=24.9 mb 
925 18.7 io 
850 14.9 d2- Ge 
700 59 4.4°C 
500 -10.7 4.9°C 
400 REO DOPE 
300 -37.3 
250 -~46.2 
200 -~56.2 
150 ~67.4 
100 -80.7 
9999* TL =22.9 


*Code 9999 indicates data taken from the top of the constant- 


flux level (CFL) of the FNWC initial data program. 
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Table I(b). 


Example of the corresponding radiative sounding 


with mixing ratio listed at odd k-levels (Fig. 2). Addition- 


ally, water vapor and CO, absorber masses, 


cloud amounts, 


CL(1) and CL(2) where applicable in the radiative sounding. 


Pressure (mb) 


1000 
900 
800 
700 
600 
500 
400 
300 
200 


100 


GeO) 


256 


17. 


eZ 


= 1M) - 
ae 


=Oie 


=o). 


=O 0. 


8 


I, 


Mixing Ratio 


e/ke 


15.49 


12.76 


6.12 


(bevel 
1.37 


~54 


02 


00 


oa 


Absorber Masses 


Water Vapor 
(gm/cm*) 


(LAs eee 
3.349 
3.633 
3.679 


3.680 


3.680 


CO2 


(cm/cm?) 


45. 


83. 


113. 


134. 


138. 


143. 


53 


53 


35 


00 


67 


35 


Cloud Amounts 
ei CL) seb C2) 


Cla(2,) 
1.000 


CL(1) 
~650 






humidity values, esi? even though the computations for oe 
made in Sections V and VI employed other parameter-values 
for these computations. 

Mimecdcieron to the parameters listed in Tab¥e I (ay~b) ,j 
surface geostrophic windspeeds were also calculated. In the 
procedure for calculation of the windspeeds at the gridpoints, 
use was made of the 100 mb FNWC analysis to determine the 
Gmmeour heights as D-values. At certain locations, climatic 
geostrophic windspeeds were used (Atkinson, 1970). A 
ieee complete explanation of the use of the geostrophic 
windspeed calculations and its use in connection with the 
planetary boundary layer sensible and latent heat transports 
memeeound 1m Section V, and the windspeed values are listed 
mieetimat section. 

forpersorm the radiative calculations indicated in Table 
jo), data-processing methods were applied to convert both 
the temperature and mixing-ratio data at k-lWevels to obtain 
eiemecloud amounts, CL (1) and CL (2) at certain required k- 
Sevels boundaries (Fig. 2). It should be observed that all 
Seumcaings in this study start at sea level with the approxi- 
mation of surface pressure T = 1000 mb. Thus the k-levels 
Semrespond approximately to the FNWC o-levels and are taken 


Sevier U0. , o0U., ....,,200., LOOT, 0.0 respectively. 


A. INTERPOLATIVE PROCESSING TO K-LEVELS IN RADIATIVE SOUNDINGS 


Phono c ta Glnes 
AGeeden intandatory Levyel of Table I(a) between the 


m0, ..., LOO mb, gridpoint temperatures were listed. An 
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assumption was made that at the top of the atmosphere, the 
temperature was isothermal from 100 mb to p = 0.0 mb. The 
listed FNWC sea-surface temperature was taken as Tio: At 

the remaining k-levels, air temperatures were derived from 
either the corresponding FNWC listed temperature or where 
Meeessary by interpolation from listed values of the original 
FNWC data. The interpolation scheme used was that of the 


mec -point Lagrangian formula in the vertical. 


T, = T, ((Py-Pz) (Py-P2))/ CP - Py) (Pg-P)) 


+ Ty ((P,-P,) (P,-P) ACP, -P,) (Py-P)) (2-1) 
Ren (P Pe) (PPC (PE Pee a- Pade 


fmetsenoted that in Equation 2-1, k is the required pressure 
tower needed for the aoe Mpternc Lat dOlemeNe sehow.e 1 Po lies 
above Py andmene spressume levels Pi and P., ive= bellows ae ved 
k. The temperatures Ty» Ty> and T, ane at tie Same. leveds 
as Pa» Py and P. Pe seeetiey Cony 

meee iumidities At K-Levels 

Moisture parameters from the original gridpoint 

soundings (Table I (a)) were converted into mixing ratios 
at each of the original sounding levels (Table I (a)). 
Then utilizing Equation 2-2 the near-surface vapor pressure 


(taken as e,) Was employed to calculate the mixing ratio at 


fee=) 10 as follows: 
Ae ane = Omyoy (e,; ,/1000) (D225 


e a 2 . ° wet | 
where tibiae °x 2 12 mb and diooo 25 17 gem(kgm) . The 
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remaining mixing ratios were calculated from the original 
FNWC dewpoint depression data levels (Table I (a)) using 


EM@itatdon 2-5 {after Fleagle and Businger, 1963): 


-B(T-T,) 
Bol. 972 A-B/T D : 
ee lr} (2-3) 
where q(P) = mixing ratio at pressure level P and 


A = 21.656 
B = 5418.0°K 
T= temperature (deg. K) at level P 
Tp = dewpoint (deg. K) at level P 
T-Ty =—"dewpoint depression at level P.- 
iemeonpuccad q-Valucs were subsequently interpolated to k- 
tevels using the procedure indicated in Equation 2-1. An 
peammole Of the interpolated q-values is shown for the case 
memene radiative sounding at point (1, 1) (Table I (b)). 
3. Humidity Extrapolation to p ¢ 300 mb 
DUCRCOUtNemtactetnat most radiosonde humidity data 
for p oOvU MDE tSwetener tnreliable or no€ available, a 
power-law extrapolation formula (Equation 2-4) similar to 
one used by Smith [1966], was employed to obtain q-values at 


Iw—o5>, 2, and 1: 


q(P) = nae A _ 
qc ( 500 aah 
merc q(P) = mixing ratio at level k in g Kg * 


dc = iene Gatto at level P —- 500 mb in g Kg 


re—moLcostncn tem cvye lai in mb, 
Equation 2-4 was solved for the best-fit parameter A, relating 
the variables y = log(q(P)/qc) id= —lowth, SOO} SIX ¢G- 
values from the original sounding including P = 1000, 925, 
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850, 700, 500, and 400 mb from Table I (a) are used in the 


Bese tit solution ({e.g., Crow et al., 1960): 


Wes 


ee (2-5) 


flo ote pate 


We i 
re OR O 


foeen the determination of A from Equation 2-5, extrapolation 
of q-values from P = 400 to P = 300, 200, and 100 mb, was 
performed utilizing Equation 2-4, 

The procedure of the previous paragraph was carried 
Siemat cach €ridpoint. As a test of the usefulness of the 
A-profile technique in determining q-values above 300 mb, a 
sommeelation coefficient was statistically computed for each 
Meemomrie. The correlation coefficient computed was between 


y and x and had the form: 


bbs pe 
Mout} ues i 


— Oe ON 


(2-6) 





pas 


As was found in the studies of Jenks [1974] and Warner [1974], 
fyewcorrelation coefficient was found generally to lie in 
Mmmemrance of .95 to .99. Such high correlations indicate 
mic ne A-protile technique for determining upper atmosphere 
q-values has realistic estimation value for the vertical-scale 
involved. 
4. Pressure-Scaled Absorber Masses 

By using the mixing-ratio values computed at the 

five odd k-levels, the layer pressure-scaled water vapor 


absorber masses denoted by AU (Fig. 1) were calculated for 


Sul 





each layer centered at the five odd k-levels. The general 
equation for computing the pressure-scaled water vapor mass 


Was. 


qj, 4P Be +72 


WORSE pley S—— | aes 


(2-7) 


where the pressure ratio raised to the exponent 0.72 has 
been employed (after Moller and Raschke, 1964). k is the 
odd-level centered in each vertical layer concerned and 
AP = 200 mb was used for each of the five layers. The pro- 
cedure was followed for each layer (8, 10), (6, 8), (4, 6), 
Weed), and (0, 2). 

The algorithm for computing the integrated water 
vapor mass above the earth's surface was then developed as 
follows: 


U(10) = 0.0 


Wie =) Vial O))) eee heat 0) 

BiGGy = — UGS) aU CG ces) 

U(4) = U(6) + U(4, 6) (2-8) 
ip eU4) = U2, 4) 

U(1) = U(2) + 0.5U(0, 2) 

Wem (Zn (0, 2) 


In the computation of the carbon-dioxide scaled 
absorber mass, a computational routine similar to the one 
used for water vapor was employed. Because it is customary 
to state the CO. aeSopPDememlasses In terms of the N. T. P. 


volume of the gas in cms over a cm’? of the earth's surface, 


the CO. absorber mass (AC) centered on an odd k-level becomes 


B32 





(2-9) 


» ONE Pk 0.72 


AC, = C(k-1,k+1) = 3.14 x 10°" ( 0 )( aes ) 


Here k is the odd-level indicator, 3.14 x 10° is the con- 


Stant mixing ratio of CO, by volume, while AP/ gp. 25) thew. 


2 
lmeer. thickness of any layer of pressure-increment AP = 200 


mb. oe is the constant density in a homogeneous atmosphere 


of depth H, in which 
P surf) = 273.16 deg. K and ne Ouse ominbe 


Famaity, Equation 2-9 1s shown to be represented in the form 


P 
k ae 
) 


=a} AP 
3.14 x 10 (H)C a973-795 )C TpTs77s 


C (k-1,k+1) 


Ce 110) 
RT 
H = — = 7) OOS se O09 ei. 
Hemlatzon 2-10 was then utilized for each odd level k = 9, 7, 
Bee, 1. 
As was done in the water-vapor case, an algorithm 


mmeear tO Equation 2-§ 15 developed to find the integrated 


G@aupen dioxide mass above the earth's surface: 


Gel — ono 
Cee re (Cone 1@@s 10) 

Ge) =Taey + C06, 8} 

GE = 0G) > Mee 15) (Zen) 
EQ Sc ew, 4 

Cee = nC) + 0.50. 2) 


Gay = Ceres 2). 





B, CLOUD PARAMETERIZATION 


The relative humidities and thus the saturation vapor 


pressure at levels k = 5 and k = 9 are used 
moms Of the £Lractional cloud cover amounts 
and (8, 9). The saturation vapor pressures 


were computed using the following equations 


Zit!) : 
621.97 exp [A- => ] 
5 
qa, (5) = 500 
CC iee xp | A = — ] 
9 


thes OE 500 


where A = 21.656 and B = 5418.0 deg. K. 


in the calcula- 
in layers (4, 6) 
at each level 


(Gia, EGiat Lon 


(2-12 (a)) 


C5 PAOD 


Mie relative humrdities were then calculated from: 


(5) 


It 
OQ [-O 
—" 
rc 
wn 


Reel. \Go:) 


q(9) 


R.H. (9) a. (9) 


(2-13 (a)) 


(GUS 101) 


Meemieetnese relative humidity values, the fractional cloud 


amounts for each layer were parameterized employing equations 


Star to Smagorinsky [1960]: 


Cin n=" OMOCRECS)) = 0.7 


CID) 


i 


Be3G MORuICON ee 2.0. 


ee) 


(2-14 (b)) 


Equations Z-14 have been employed in the radiative studies of 


Jenks [1974] for a spring day in 1973 and by Warner [1974] for 


mmease in mid-October, 1973. In the current study, Equations 


Meer eore tested in a winter sample of FNWC data for 16 


January 1974. 
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Smagorinsky's parameterization of CL(1) and CL(2) per- 
mumes cloudiness fractions of 1.30 and 1.333 respectively with 
RH = 1.0. Values of CL of this magnitude were allowed orig- 
inally by Smagorinsky to suggest the type of supersaturation 
which accompanies precipitation. However, in all of the 
radiative models in use (that of FNWC, UCLA, NCAR, etc.), 
mwas Considered correct procedure to limit CL < 1.0 in 
both layers. In the radiative calculations of this study an 
arbitrary reduction by one third of that computed by the 
Smagorinsky formulations (2-14) was made so that neither 
fractional cloud-cover could exceed unity. For radiative 
calculations, the degree of supersaturation in clouds was 
eonciadered not of importance. 

Thus CL(1) and CL(2) computed by (2-14a, b) are 
mpeeaced respectively by their comparison cloud-covers, 


were)” and CL(2Z)', defined by 


(es 27 5) GINGE) 
(2-15) 


CEG 


Zo, Glee ) 


mebed at Cach gridpoint. The objective for testing the 
meauced cloud-covers (2-15) was to determine the relative 
eeimce Of the displacement of the radiative balance of the 
earth atmosphere system. Such radiative tests were esvecially 
crucial since evidence from Warner's [1974] tests made with 
Mamoditied (CL(1), CL(2)) indicated that the planetary albedo 
was larger than observations from NIMBUS III (Raschke, et al. 
M75). the one-third reduction in the Smagorinsky calculated 
amounts used in this study is an initial effort to tune cloud 


amounts to give radiational results in closer agreement with 


3D 





mene latese satellite*results such as are found in Raschke 
eeeat. (1973). Aliso, the fractional cloud amounts were 
considered functions of large-scale effects only. Small 
scale convective activity, seasonal conditions (other than 
the soundings themselves), and latitudinal effects were not 
considered in specifying this reduction factor. The purpose 
of the reduced cloud-parameterization modification was for 


estimating radiational effects only. 


Gee wINDSPEED USED IN TURBULENT HEAT TRANSFERS 

For all gridpoints located between 5°N and 5°S, gradient 
windspeeds were taken from the climatological data for the 
month of January (Atkinson, 1970). Gradient isotach values 
Menemutilized at the tropical gridpoints noted due to the 
email Coriolis parameter £ in Equation (2-16)(below) in the 
latitudes between 5°N and 5°S. 

prea other erildpomnts, Equation (2-16) was utilized to 


Compute geostrophic windspeed given by 


pe ee g lazl CP aalis 
BFA 
m 
here g = 9.80665m sec 
= De sary 
m = map factor = T + Sing 


latitude 


“- 
i 


d = 381. x 10?m = nominal FNWC grid-dimension 


4 ; 
( 36-400 ) Sind (2-17) 


| Az | oan (25 - 24)- (2-18) 





ie z-ro) wand (2-18), Az 1s the Contour-difference cen- 
fered in general on the gridpoint. At the map edge grid- 
points of meridians l, 1, ee fae was possible to compute 
only a forward differenced version of Equation 2-18 (cf., 
eee oo). Thus at edge-points "0", Equation 2-19 was used 
(below). Note, the grid distance d/m in Equation 2-19 is 


one-half of that used in Equation (2-16). 


Mel = ghz oils) 


with 


Az = Due Zee (25 2825) 2 « (2-20) 


Pueeececostrophic windspeeds were then converted to surface 
Windspeeds using the following empiricism due to Langlois 


and Kwok, [1969]. 
V = 8 Vi 2.2 (2-21) 


where ue = surface windspeed (mps) 


- = 1000 mb geostrophic windspeed. 


Meee that in the interval 5 to 10 degrees of latitude in 

either hemisphere, the Vamienor the Coriolis parameter was 
feerarily set at F = .25 x 10 “sec, following Kaitala 
[1974], while for [¢] > 10°, the actual value of f for this 


particular latitude was used. 


Sil 
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Figure 3. Section of FNWC polar stereographic meridian il- 
miovrating the method of obtaining contour gradients in the 
vicinity of a gridpoint. The values of contour-heights 
correspond to the 1000 mb D-values in the neighborhood of 
the gridpoint. 
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pee THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas for computing Emissivities for ter- 
Meseuial tlux calcullations in the atmosphere were developed 
by Sasamori [1968] in association with the operation of the 
Neen General Circulation Model. These empirical formulas 
for water vapor and CO. emissivities were matched to the 
theoretical values used in the Radiation Chart developed by 
Yamamoto [1952]. ‘This chart was proved by Jenks [1974] to 
Momteeniaadte cnoueh £0 simulate the empirical methods of 
Sasmorl and 1s used here as a schematic guide for the inte- 
gration of the radiative transfers through the various layers 
tmene radiative soundings described in Section II. 

LomeusceliMmerncmENNGemecating patkage, the essential Vong- 


wave flux formulas required are the following: 


Fig” = het k= flux at -eartn, « — 10 

Fe = Nemeth fluxes leveliek = 6 

F,* = net IR flux at level, k = 2 
Folg = met IR flux divergence in the layer (6, 10) 


F26 = net IR flux divergence in the layer (2, 6). 
MmOordcrmeeomcc ENEnC lhe thUe-aivenmacnces (1,.¢€,, cooling 
fates) in the layers (6, 10) and (2, 6), the differences 
Aa = Fi9* and F,* - a (MicmDcmcolpurcd ss lhe detailed 
scheme for making such computations for various combinations 
Seecloud cover Ch{i), CL{2) is similar to that used by Warner 


[1974] and is explained in the following subsections B, C, D, 
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and E along with the diagrams Figures (4a,b,c,d), (5a,b,¢,4d) 
and (6a,b,c,d). A complete listing of net IR fluxes and 
cloud amounts was computed at each gridpoint. An example 
Geecthrs listing is shown for "“full-cloud" amounts, Equation 
fore, in Table Ii. One must further have a physically sound 
representation of the emissivity Ec as a PUNGCELONeer berth 
water-vapor and CO, absorber masses along the sounding, that 
fPmealong the path of imtegration in order to make IR net- 
flux calculations. For further explanation of the emissivity 


formulas (after Sasamori 1968) used in the quadrature scheme, 


meocor to Appendix A. 


Bem! FLUX Fig” 

As depicted in Table I (b), radiative soundings have been 
computed as the combination of parameters U(k, 10), Clk ..10)5 
Ths imeoeewcach tcaquired level k. Also, listed at each grid- 
emi are the cloud-cover parameters CL(1) and CL(2) com- 
pled from (2-14). The grid area may then be thought of 
being composed of areal fractions (or weights) W(0, 0), 
foe t), W(l, 0), and W(0, 1) which are defined by Equations 
So |) ee wocaeas follows: 

(a) totally clear skies, case (0, 0) 

COC 0) amet 2eniC1)) Gl Ch (2)5 oa) 

(b) oOVewCcCast in both Layers, case (1, 1) 

ie LS CGE) (CER }) (3-2) 


(c) overcast in upper layer only, (case l, 0) 


Pleven ( Clot) ) (1 el (2) ) (3-3) 
(d) overcast in lower layer only, (case 0, 1) 
GO een CLC ty) (Cli(2)). (3-4) 
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By multiplying the appropriate weighting factors by 
fme respective net flux computations Fig’ ( Oa O18 Waar er 
Fi9* (0; Ij the Gomposite net £luxes at level k = 10 for 
the combined cloud (CL(1), CL(2})} coverage over the grid 
area are then obtained. Such composite net flux values 


have the following form using level k = 10 as an example: 


* = * 
ar IGM rac ON Wigan) is ats@leeed ) ~~ 


> Cle Oe ls Ole IO, LAR Os We 


ieee Net Flux with Clear Skies, Fig (0, 0) 
In Fig. 4 (a) the hatched area represents downward 


MORO ret nein eee ra le to xm: 
¥10 


oe dB ; By = oT, 


mx FI0OY at k 


\ 


Fue 
B=0 


By imine otephen-Boltzman blackbody flux at Tigo steeei 1s the 
combined water-vapor and CO, emissivity along the sounding, 

eons related to absorber-masses as detailed in Appendix A. 
Mmeemmmnatched area therefore depicts the net flux at k = 10 

With clear skies given in the equation 


Big 


* - = = * 4 “oe 
i oo sud ap (1 ed Beene So*Ts (3-6) 


ee Net Flux with Overcast Clouds in Both Layers, Fi (1,0 


rmeoUremo Dy ti lustrates this case. Byeusing the 
trapezoidal summation techniques, F10+ at k = 10 is depicted 
wiiee again by the hatched area representative of the case of 
PumovercGast With base at k = 9. The unhatched area depicts 


the integral 


15 





10 
F 


ae 


fash) = (1 -.e _)dB Gap 


10° | seh 


B=0 
with overcast clouds at both upper and lower levels (Fig. 4 


(b)). 


3. Net Flux with Overcast Clouds in Upper Layer Only, 
Fip CG, 9) 0 


Here an Fig. 4 (c), the hatched area depicts the 
total F104 emitted down from the cloud base at k = 6, plus 
aieinerement in the upper right corner of Fig. 4 (c), which 
represents the contribution to FY from the water vapor and 
CO. i the atmospheric layer from k = 6 to k = 10. The un- 
Mmemencad area depicts the integral 


Boo 


Fig, 0) = nA (cll Eo) dB. (S00) 
4. Net Flux with Overcast in Lower Layer Only F960: iL) 
This case is the same as depicted for case Fig GC: IL} 
ieeee 4 (d) and is explained in Section III.(B).2. Ser is 
Pie temperature-dependent emissivity for the completion of 
the integration from level k = 1 to the origin (B = 0, € = 
mavyeot the Yamamoto Radiation Chart, applied for T < 210 
gem, K. Ce Powone COMperatume- Independent —eMmassivity applied 
mere £ > 210 deg. K in the level from k = 1 to k = 10. For 
the specific formulations of ee and oe after Sasamori [1968] 
see Appendix A. 


fe ihe Composite Figs (CHC bmcheZ  IeEquation 
The principle involved in compositing net fluxes has 
iercaday been expressed by Equation 3-5. The following equa- 


tion can be thought of as the general form of the quadrature 
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scheme used to evaluate the integrals for computing Fi” eae 
mmaconommarcronvot (CL(1), CL(2)):. It is also used in com- 
puting the integral for the composite cloud-cover case Fig*> 


Gemresponding to (3-5): 
Fi9*(CL,,CL,) = [1-CL(2)}{(B,,-Be)-. Sle, (8,10) (B, 9 -B) 


+ (€,. (8,10) +e,,.(6,10)) (Bg-Bg) I} 


+ 


OGL )) G SHED) Bs Gey Os 0) 
2S 0D, UO) C3. aE Gy, nO} (3-9) 


+ €¢(2,10)) (B,-B,) +(e, (2,10) 


+ 


2 EO) EG SS Op pe 1 A 


+ 


CL (2) Cs pele) ldlee Sa ls Oe 


tyeesctcie nt) = Cih(z) = 080, Equation 3-9 reduces 

fomeme form needed in calculating Fig (9, 0)... Simi laaiy wore 
Can obtain the equation for computing Fi) (1, 1) as depicted 
Smpene Chart in Fig. 4 (b) by setting CL(1) = CL(2) = 1.0. 
Batally, the other Boe Stibeaee = Fig GC, 0) and Fig* (0; ig) 
eaeoe formulated by setting CL(1) = 1.0, CL(2) = 0.0, and 
mime) = 0.0, CL(2) = 1.0, respectively in Equation 3-9. 
maiiebar procedures were followed for computing net [IR flux 
mere Vels kK = 6 and k = 2, as described in the Subsections 
mete (ries. 5 (a), (b), (c), (d} and 6 (a), (b), (c), (d), 


mespectively. 


47 





me NET FLUX Fe 


1. The Composite HAG Equation 
The following expression, which represents the Be 
composite net flux, may be thought of as the weighted sum of 
the unhatched areas of Figs. 5 (a), (b), (c), and (d), with 
Ene appropriate weight factors W(0, 0), W({1, 1), WC(1, 0), 
mee, 1) of Equations 3-1, 2, 3, and 4 applied. The re- 


sulting expression for Fe*(CL,; CL.) becomes (after Warner, 


1974) 
i = [1-CL(1)] {Bp-.S{e, (6,8) Bg-By)*e, (4,6) (Be-B,) 
See onc (210) (By 8) Ce oo ieee lc Pe oee) 
een) at se l= Cl Cl Go iC teGie(2) it (By eel. 
-5(e, (6,8) +e, (6,10) )]}+CL(1) {(B,-B,)* (3-10) 
ree (058) 1) 2 CL (11) (1-1 (2) ) 4 (BE, 1-8) 
Peete go 8) = 0,1 0) T- 


2. Net Flux With Clear eee Fe* (0, 0) 
As was described in Section III.B.5 concerning the 

miewor the composite case (Equation 3-9), Equation 3-10 can 
be reduced to give expressions for ae for the various cloud- 
mmeeercases (0, 0), (1, 1), (1, 0), and (0, 1). The resulting 
expressions are depicted as the unhatched areas in Figs. 5 (a), 
mee (c), and (d). For the clear sky case of Fee, Equation 
wet0 would be employed with CL(1) = CL(2) = 0.0. This is 


depicted by the unhatched area in Fig. 5 (a). 
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See inth Overcast Clouds in Both Layers ee 1) 
For this case, Equation 3-10 would be used with CL(1) = 


mraz) = 1.0: This is depicted by the unhatched area in Fig. 
oe (b). 


feeeNet Flux With Overcast Clouds in Upper Layer Only 
F Ge 0) 


“6 
mie expression needed in this case, and which is 

depicted by the unhatched area in Fig. 5 (c), is formulated 

igeeusing Equation 3-10 and setting CL(1) = 1.0 and CL(2) = 


G0. 


S. Net Flux With Overcast Clouds in Lower Layer Onl 
Eb Onan 
“6 ) 


Micmac NOneMce aed sine tilts) Case wmrepresch tame), 
the unhatched area in Fig. 5 (d), is formulated by using 


Pomation 5-10 and setting CL(1) = 0.0 and CL(2) = 1.0. 


pe NET FLUX F,* 
A procedure analogous to the ones used in the computa- 


j@pons of F.,* and ae Mec fluxes ean be Used in Calculatiren 


1 
meme t flux at level k = 2. F,* ismerepresented pictorial ly 
by the weighted sum of the unhatched areas of Figs. 6 (a), 

(b), (c), and (d). Once these areas have been computed with 
miieemroper weighting factors applied, as was done similarly 


mmeecduations 3-9 and 3-10, the following general expression 


for the composite result for F,* es wlligs 


Po = Peeieh GUE ye ae (24) (By Boe (2,4) 


+ €.(2,6)) (Be-By)*(€,.(2,6) te, (2,8))* 


(By Bere. (1,2)(B,-B,)te__((0,2),T,)*B,)1} (3-11) 
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eleingt) l-CL(2) )*1 (35 ebey lien (251) (3-11) 
+ €. .(2,10))]}+CL(1) (B,-.Sfe, (2,4) (By-B,) 


+ € (1,2)(B,-B,) +E, ((0,2).7,)*B,]} 


Poametivon 3-11 can also be derived in the computation of the 
meer luxes for the general cloud-cover cases once the proper 
weight factors for the simplified net fluxes are introduced. 
These simplified F,* Gases ecomLeoponanstoetcne (0, 0), (1, 1); 
feeeeo), (0, 1) cloud-overcast/clear combinations and are 
mepresented by the notation F,* (0, Oy Fo*(1, sh a Es Oy 
F,*(0, Meerespeemively, Whe correspondme chatrtsdepmerlons 


paemalsplayed in Fig. 6 (a), (b), (€¢c), (d), respectively. 


Pome PLICATION TO HEAT BALANCE COMPUTATIONS 

inmerder tomcompute the heat balance of the atmosphere, 
layers (2, 6) and (6, 10) are treated as subject to solar 
feaerie (Section IV). In order to complete the atmospheric 
radiational balance, the long-wave mean cooling rate in the 


FNWC model was computed as follows (after Kaitala [1974]): 
F26 = F,* - FL* (5-12 (aD 


et 0 


I 
al 

% 
zo 


6 10° es FZ (oe) 
Wieme the combination Z6 and 610 indicates the layer boundaries. 
mate1ons 3-12(a) and (b) represent the cooling effects by ter- 
Pestrial radiation in the indicated layers, which is applicable 
fmeine time of the analysis. The radiative balance at the sur- 
Bace is computed using the absorbed solar insolation (Section 

ie 5-C) and the net flux Fig* PmOmenalattons(5-9). /A conplete 


heat balance at the surface is computed by subtraction of the 
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evaporative and sensible heat transport (Section V) from the 
miecace radiative “balance.'’ All these balances are dis- 
played in cross sectional display forms in Section VI. The 
IR-cooling analysis was carried out at each gridpoint and 
the cooling rates were stored by gridpoint along each meridian. 
The zonally-averaged atmospheric balances are detailed in 


Seecions Viland VIII for 16 January 1974. 


F, STATISTICAL RESULTS AND COMPARISONS 


A statistical test of the F,,* numerical results (Equa- 


10 
meron o-9) was conducted by utilizing a linear regression 
Program tor regressing the 93 values of Flo’ Gnetnien con 


responding sample values of B and Bigve aise predret ors me dine 


10 
program used was based on the BMDO2R in the Biomedical set 


of programs (after DIXON, 1973). The resulting form was 


Fin* = Big (a + dre) Gauls) 


where Bio = Stephen-Boltzman blackbody radiation at the sur- 
face with temperature Tip: Bio = Sta Tio! 


surface layer vapor pressure (mb). 


HI 


e 
meme > 1n Equation 3-15 are the desired regression coeffi- 
cients and are dependent on B,,) and Bi gVve (aber Brune 9 52). 
mies predictand used for the regression was Fig’ Wlgiaee}ol aug 
cluded the downward contribution of CO, as well as that of 
fewer Vapor in the computation of net flux at level 10. The 
Seaeistical equation found in this study is Equation 3-14 


where Ra PE eclomimulea LecmeCorreration coetrici1ent . 
* 
Fi0 


R 
m 


S8GL8.. 5 O20, BE 
10 10 (3-14) 


oO Soi 


aya 





19 and By yve are good 


Bemeedictors of Fi” and therefore of Fa obtained from 


The high value of the ae Shows that B 


— ; ae ; , 
Foo Bio Fa: Pie ROC Gt dO aon: Foo HMO os pean Ss 
feutcion of the latter for the downward flux at the surface. 
F is given by 


F 


q Pela apres 


or (3215) 


F Bi, (. 6436 NAS 


d 
Pomacton 35-15 1s to be compared with a similar result after 


Swinbank [1963] 


aa = Bop lnoe + ,037Ve). (3-16) 


memecoummary, the high value of the correlation coefficient of 
Paieation 3-14 along with a strong similarity between Equations 
3-15 and 3-16 seems to indicate that Equation 3-14 may be 
regarded as an oceanic version of the Brunt downward- flux 
ealaeron, Also, it is noted that i depends most strongly 


Hoon B,, and only slightly on (Bi Ve)-. This was to be ex- 


10 
meeted over the oceanic regime where nearly constant relative 
Mumadities exist over large, expansive areas. 

Pmescecond Statistical test was based upon a possible re- 
jationship between cloudy sky and clear sky cases of Fig’: 
Imeethis test, the ratio Fig (CL)/Fi9* (0, 0) was used as the 


Meodictand. The predictor was an approximation to the com- 


Peotte total opaque cloud at each gridpoint, defined as: 
Cl eee) tern 2) (le > GI C1) i. CS) 


Mie regression was run twice, first utilizing Fig” (CL) com- 


pered using the amounts of CL(1) and CL(2) as calculated from 





the Smagorinsky formulas (Equations 2-14(a), (b)) and second, 
tea 11 Z2ing Fig” (CL) computed from the two-thirds amounts of 
CL(1) and CL(2). The resulting regression formulas based 
upon a best-fit to the 93 soundings used in this study were 


meeoced in (3-18, 3-19) 


ech hy (clye-eE C0, 0) (1-.75CL},eR =—n9916 (5-18) 


Benes (CL) er D, OVGle, FSC Se GO Gieils) 


* 
10 ( 
Pmeomparison equation after Sverdrup [1942] is included as 
(3-20) 


PaqcGl) 2% O- OMG: Gel (=200 


10° § 
Gmemnieh correlation coefficients found with Equation 3-18 and 
Seeeoealong with the similarity to the Equation 3-20 (after 
moemeanin, 1942), indicates the general capability of the IR- 
madtative model employed in this study to account for the 
paeeetive net radiation at the surface, Fig (CL,,CL,); both 
mmiene completely clear and cloudy-sky cases. 


The mean values of the parameters in Equations 3-18 and 


Selewover the 95 gridpoint soundings were as follows: 


Mie GLOUD STATISTICS REDUCED-CLOUD STATISTICS 
Fyp*(CL) = .0749 Ry min” Eyre = . 0950) 2y min 

7 2 ee oe : . =] 
Fo im 0) = 2~WSheery min Fy" (0, 0) peomleGee y aatden 
CL = .6951 CL = .5190 


mime above table, superior-bar symbols denote sample means. 
These results indicate that the surface net flux is decreased 


by the respective ratios (0709/25 1S) = 507 with CL = 


.6951 and (.0568/.1518) = .3740 with CL = .5190. These results 
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midweate tmat the downward flux F,,+(CL) 1S, on the average, 


10 
rey  peTecenlt larger than the clear-sky case utilizing the 
eloud-parameterization model (Equations 2-14(a), (b)) in 

this study with full cloud amounts. On the other Fig? (CL) 


moe the two-thirds cloud-model, the downward flux is increased 


Bempenly 37.4 percent relative to the clear-sky case. 


Ga COMPARISONS WITH CEIMATOLOGICAL RESULTS OF F, 


Comparison was made of F, by the model-computations with 


2 
satellite measurements of total long-wave flux to space after 
Raschke et al. [1973] for the NIMBUS III period 21 January - 
3 February 1970. F, dseused instead of F,* Since the satel- 
Mmeeenaoes Not measure the downward flux at k = 2. F,is 
computed using Equation 3-11, which was also used for F,* 
burewomitting those emissivity-terms involving levels above 
k = 2, that is Ec (42> IL have Ec (2, Oe enone Lala t lol eo —lele 
Miteec Lerms are associated with downward flux at k = 2. The 
meoees tor 16 January 1974 were computed both for the full- 
cloud and the two-thirds cloud-amount cases, based upon all 
Sr1dpoint soundings. These EF - values were then interpolated 
memnnole multiple of S° latitude along each of the four merid- 
mages and the results latitudinally averaged to get a "zonally- 
averaged" value valid for each 5° latitude band used in this 
Seudy . 

Table III (a) shows the zonally-averaged F, - values com- 


bared with those after Raschke et az. For consistency in 


comparison of the F, results deduced here with those reported 


Hi 
by Raschke, the latter's results were obtained by interpola- 


tion from charts using the same oceanic meridians as those 
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Radiative Statistics | Cloud Statistics 


F2 Fo Fo FULL-Cloud 2/3-Cloud 

Lat (FULL-CL) (2/3-CL) (Raschke) CL AVG. CL AVG. 
205 Pros laa ro 2o m3 0 BOS ~758 
ILS) a .334 mee ~927 Perio 
10 . 344 . 354 F OW it .867 .634 

a) epee: 5 OOD ~ 347 . 746 D0 

0 vo6l . 369 cioul . 790 .5596 

= eo tal oS 5 SeNe: . 766 oe 
10 = Sie Aye iiss 7/ . 388 Omg OO 
15 S3O0 Zoos . 394 TOO 2 ~342 
20 BoOU . 388 paigak .443 P24 
25 mei ei OOZ Soil .463 ele, 
30 .354 eas: . 344 009 ~411 
SUS) 7336 . 342 .319 POG Wood 
40 Polke 2S IS, . 304 .628 ESV, 
45 2275 PAL) - 292 TOO . 740 
50 . 264 (PRES - 210 ~911 5 eS 
55 PAR e253 i Zao ~931 . 749 
60 . 209 . 226 AS . 764 “Be 
65N 200 S200 20 Oz 0.0 
ae 333 . 338 344 688 514 
Avg. ; 
Table III(a). Comparison of zonally-averaged longwave fiux 


memepace, F,, aS found by this study for both full-CL and 
2/3-CL cloud cases for 16 January 1974, and by Raschke et al. 
(1973) based upon NIMBUS III measurements. Also included are 
composite cloud amount fractions used in the two cloud 


models of this study. Flux values in ly min7?. 


B19) 





used in this study. For further insight into the compari- 
moms, thie cloud statistics, CL and 2/3-CL values are listed 
at each latitude, where CL is given by Equation 3-17. In 
wae bottom line of each column in Table III (a) is listed 
the cosine-weighted mean of each set of the column values 
over the listed latitude range. The resulting oceanic-means 


one EF, Speiimea Dyetne Varlous techmiques are as Lollows: 


F, (Full cloud amount this study) = .333 &y min 


F,(2/3 Cloud amount this study) Oo Oomy: min ” 


.344 2y min 


F, (after Raschke) 


Nete that the Raschke results do not correspond to speci- 
je known CL - values. The only conclusion that can be drawn 
iogom table IlI (a) is that both cloud-amount systems give F 
meeults which are reasonably close to the values reported by 
masenke. This is especially true in the subtropics, 15°-20°N, 
where peak values of F, are found, and in high latitudes where 
minimum values of F, (after Raschke, 1973), are found. By 
these results alone, neither the full-cloud model nor the 
2/3 CL can be selected as preferable. The limitations of 
the comparisons made here are obvious, when it is realized 
that at latitudes between 20°S-5°S and between 60°N-65°N in 
momre Ili (a), fewer than 4 meridional lines enter into the 
zonal averages listed. For all other latitudinal average 
Values, four meridional lines were used in this averaging. 
m@aese same limitations will apply to Table III (b) below. 

Note finally that the model results F, (full) and PF, (27 5 aCe) 


mee based upon a single-day's observations, 16 January 1974. 
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Pao Oo Gas By ot 





Lat CHU Ey Cir) (27 oe el (Budyko ) 
205 .054 .078 .O91 
15 .050 7 O79 Ser 
10 .055 . 086 ake 
+9) 7063 .090 .085 
0 .052 .082 .084 
) .056 . 086 .080 
10 .085 . 107 . 086 
15 .094 LES .094 
20 - 102 ~120 LOZ 
295 .104 ~120 .098 
30 . 096 Pri Ws .099 
35 .107 PLZZ ~ 102 
40 . 096 -108 . LO9 
45 .056 . 067 .098 
50 .054 . 069 099 
OO .052 .068 .102 
60 .085 OG .104 
65N 2095 -095 .095 
Wt. Avg. .075 .096 .094 
Table III(b). Comparison of zonally-averaged net fluxes at 


the surface F,,)9* as found by this study for both full-CL and 
2/3-CL cloud cases for 16 January 1974, and comparative 
Values aiter annual climatology of Budyko (1956). Flux 


values in ly min™’. 


Sil 





H. COMPARISONS WITH CLIMATOLOGICAL RESULTS OF Fig” 

Table III (b) depicts the zonally-averaged values of 
cla Fig” at the earth's surface obtained in this study using 
both full-cloud model and the two-thirds cloud model. A 
Somparason of the model-computed Fig*-values Hoe telvey awl 
and two-thirds cloud cases is then made with annual mean 


values of F,,* obtained from climatological studies (Budyko, 


10 
iso). The Fig” values shown in Table III (b) were computed 
meen Equation 3-9 using the full CL.,CL,. values or the re- 


aaa” 

duced (2/3CL,, 2/3CL4) im the respective cloud-models. 
pamec) bUdyko's mean values of Fio® were based on annual cli- 
matology, his overall latitudinal range of Fig” Viaies as 
not as great as that of Fig” by either cloud-model presented 
mielaple DIT (b) for 16 January 1974. 

It is to be noted that the locations of highest values 
of Fig” after Budyko agree with those computed for both Fig* 
eel} and Fig (2/3 Cle namely iit themsubtuop1es (latitude 
20°N), and again in the latitude range 35-40°N. The "hemi- 
Spheric"’ weighted mean of ae = OSS ay min One eine es 
duced-cloud model is slightly greater than the cosine-weighted 
mean of Budyko (.094 xy min’). However, the weighted mean 
oe Fi9* (CL) resulting from the full-cloud model is consider- 
ably smaller, Fi Gri =. 07S. 2ky min’. The chief differences 
in Fi 9% (CL) as a function of latitude compared with that of 
10° Cor- 
responding to the relatively high composite CL-amounts in 


Buavyko 1S accounted for by the marked minimum in F 


Mow-latitudes (from 20S to SN) and again in the latitude band 


45-S55°N where in both areas typical minimum values of Fy * (CL) 


58 





= 055 Ly min. occur. Such a minimum of Pig” Pile ie wc eairOns 
45-55°N is generally associated with the location of the 
Boelar-front zone. However, in this same region the 2/3-CL 
model gives average values of Pio” = .068 Ly min’. The 
Zomal-mean values of CL and 2/3-CL applicable to Table III 
(b) by latitude are exactly as listed in Table III (a). It 
is not clear from the comparisons of Table III (b) whether 
the higher values of Fi9* deduced from the 2/3-CL model 
represents a closer approximation to reality than the Pig” (CL) 
fees resulting from Equation 3-9 with the full-cloud model 
Walues. 

It is to be expected that analogous comparisons made 
using the solar-insolation model (summarized in Table IV) 
imolietindicate some relative advantage for using the full or 


two-thirds cloud model. 





IV. SOLAR RADIATION - 


A. PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k = 0 
(top of atmosphere) was 2.00 hy min (Joseph, 1971). This 
meee Value was further depleted by 4% to account for the 
attenuation caused by oxygen and ozone above the tropopause. 
This left the value S = 1.92 £y min. at level k = 2 to be 
used in this study as the effective solar constant. 


Bouation (4-1) was then utilized to compute the effective 


solar insolation at the tropopause (k = 2) as follows: 
T ee 
PO = Ss ] a ] Cos Zz (4-1) 
m 
where S = effective solar constant at k = 2 
Gos z = cosine of the zenith angle for the Julian date used 


yA = ratio of the actual earth-sun distance to the mean 
earth-sun distance for the Julian date used in 
L)oaL GS, Gyrgibuehie: 
miewomtthsonian Meteorological Tables (List, 1958) gives the 
ratio yt and the solar declination 6 for 16 January, OOOOGMT, 


meeeegese values are employed. These values are, respectively: 


clio SFIS 


== 
Yr 
m 


6 -21.08 degrees of lat. 


®)i1s used in evaluating the cosine of the solar zenith angle, 


meiven by 
COsaze-  Sifieorotn « + Cos Oo Cos 6 Cos h (4-2) 


where » is the latitude and h is the hour angle of the sun 
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molative to the meridional data lines. For example, Fig. 1 


makes it clear that for 


lene. £ h = 55° 
line 2 le = LE 
lines 3, 4 | ie 


metne times of the synoptic charts (QO00GMT and 1200GMT, 
respectively). Sin ¢ was computed according to the standard 
Pemar stereographic projection formula applicable to the 


FNWC base chart given by 


i (SS [ion | 
—— ———— (4-3) 
Ga (LESAN (elas) =) 


where Tro = 973.752. Thus Sin ¢ assumes the following func- 
tional form in terms of the FNWC grid-coordinate I (Fig. 1): 


mines 1, 3, 4 


ie USES 1) ; 
973.752+2(S2-1)2 (od. (Ca) 


Sin 9 
Line 2 


O78, 782 (B21 e ; 
ee oa So (7 (22) 


or, conversely I is given in terms of 9 by 


; 2 : Cos ¢ z 
Lines 1, 3, 4 I = 32 - 22.065(;—° 2] (4-4 (c)) 
Line 2 2 82 2 Sh 205 (ee eas 
a eC raa 


li. , o2o0t Omebi nes: 1,5. -2 
I = (Ree. wiz outor Line <5 
Cerio So nO alate 4. 


Meme lines £, 2, and 3 the gridpoint soundings correspond to 


QOOOGMT 16 January when solar noon occurs at the iso0°8 
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meridian. Line 4 gridpoint soundings correspond to 1200 GMT, 
16 January, when solar noon was at the Greenwich meridian. 

A very simple partition of solar insolation was utilized 
mimcthis study after Joseph [1971]. It consisted of dividing 
mer insolation F(Z) into two parts at level k = 2. One part 
was considered to include all wavelengths A > .9 um where 
absorption by water vapor and carbon dioxide bands are the 
Mest prevalent attenuation processes in clear air. This part 
of solar insolation was listed as F(A) energy and considered 
Gueject tO water-vapor absorption but not to Rayleigh scat- 
Being. For those wavelengths A < .9 um, absorption of the 
solar insolation energy by water vapor was considered neg- 
ligible. This part of the solar insolation band was listed 


as F(S) and was subject only to Rayleigh scattering attenua- 


muenmeim Clear air. the two partitions are formulated as 
follows: 
F(A) = .349 F(Z) (me a) 
pS Sy) ie To) fs) MI BC) (abe ss Gen) 


Mienis study, the introduction of two cloud decks produced 
meaua-reflectivity effects upon both the F{A) and F(S) solar 
energy insolations. However, in the clear areas around any 
memapoint the absorption-attenuation only applies to the 
F(A) insolation while the Rayleigh scattering-attenuation 


mepiies to the F(S) insolation only. 


bee DISPOSITION OF F(S) INSOLATION 
Ine thewtLreaiment o£ the P(S) ansolation, Joseph [1971 } 


Hetermined that Rayleigh scattering reflectance in clear skies 
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/feeter Coulson, 1959) could be effectively approximated by 
least squares in the following form 
a(R) = .085 + .25074 [log ( a Sec z )] (4-6) 
O 
where Pe Seco Own Equation 4.06 T/P = 1 in view of 
mmemtact that mean sea level pressure 7 is close to 1000 mb. 
Also 


Sec z = (Cos z)° 


Tee COS Z Piven by Equation 4-2. 

The surface aibedo a(G) is another reflective parameter 
utilized in this study. Over oceanic areas the following 
mommuka tor a(G) after Gates et al. [1971], was utilized: 

ee mMaxel.06% .060" .54.(a7 -iCosiez)}: (4-7) 
ieeclear Sky Case | 
In the clear sky (0, 0) case the F(S) insolation was 
muegected Only to Rayleigh scattering reflectance a(R) and 
to surface reflectance a(G). Considering the likelihood of 
@ succession of multiple reflections between earth and 
memesphere, the F(S) insolation actually penetrating the 


eee’ Ss Surface after scattering 1S given by 


mona, ON = FUS) Poot RN 1 tlteqRe (Gite. (o(R)e(G))_ 


eee (G)) (4-8 (a)) 
Minat 1S, by 


ISi0(0,0) = F(S)[1-a(R)]}[1l-a(G)]/fl-a(R)a(G) (4-8 (b)) 


2 Cloudy -Sky Cases 


In the three cases in which clouds were present, F(S) 


insolation absorbed by the ground at each gridpoint was 
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Computed using a variation of the following equation (after 


Arakawa, 1972): 
ape ee Gal 2 (2 Gio (G) ) 


AS | ORIG) 81 2)) ROCA CG) ACD es Ce) 


ANG elei!  ~ ao 


Pema cated by the notation (1, 1), denoting CL{1) = CL(2) = 
ieee Equation 4-9 is the formula used in calculating F(S) 
mescolation absorbed by the ground in the case where overcast 
memomas are present at both levels of Fig. 2. Also in Equa- 
eon 4-9, ‘constant cloud-reflectance values were chosen, 
Mamely R{1) = .54 for the mid-level clouds between k = 4 

and 6, and R(2) = .66 for the low-level clouds between k = 
Pam. Both cloud-reflectance values are as suggested by 
C. D. Rodgers [1967]. 

For the other cloud cases, the following changes 
maweemepplied to Equation 4-9. In the (1, 0) case (CL(1) = 
mater (2) = 0.0), consider that R(2) = 0 in (4-9), from 
which it follows that | 

OC Ome el Gy) Cl aaGGy) 7 [tek (1) o(G) 12 
Mmimeene case (0, 1), R(1) = 0.0 in (4-9) so that ala 

Po CO ly) = PCS USCA ESOC Ye Lt CON SOT eas 

me Composite F(S) Insolation 

Pamaieronsm toi 4-9) (4-10) and (4-11) were 
Mamlrzed in the computation of the clloud-weigphted F(S) in- 
meeacion penetrating the earth's surface considering felare 


pecal-weiphts of the cloud combinations denoted by (0, 0), 


ae) eo ean (Oe I} eabout a gridpoint. The resultant 
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me insolation penetrating the earth's surface denoted by 


mod is therefore expressible as 


MCC ame Cte) = lodcnon 0) Wane an 
+ IS10(1, 1) W(1, 1) 
(a= 12) 
+ IS10(1, 0) W(1, 0) 


+ IS10(0, 1) W(0O, 1). 


Here the weighting factors W(0, 0), W(1, 1), W(1, 0) and 
feet) are as computed in Equations 3-1, 3-2, 3-3, and 3-4 
meopectively. Note finally that the part of F(S) insolation 
meelected to space is found by subtracting ISI0(CL(1), CL(2)) 


inom F(S). 


Pee loPOSITION OF F(A) INSOLATION 
vem ractionalsportionsot the solar insolation subject 
to absorption by atmospheric water-vapor and carbon dioxide 
mucmeovered in the following subsections. 
1. The Clear-Sky Case (0, 0) 
The Manabe-Moller absorptivity function provided the 
necessary absorptivity values for the key layers in this case. 


Mieertorm of this absorptivity function is 

Ceo) | ieee )Scemeae (4-13) 
Hemwena 1S the absorptivity applied to the pressure-scailed 
Maer Vapor mass between levels 2 and k (Fig. 2) along the 
Meee Slant-path angle z. The resultant absorbed insola- 


mmomal energy*in the particular layer (2, k) is then given 


by the Manabe-Moller relation 


Ne Moe nOec Irie sec zl. (4-14) 
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The two layers of interest in which absorption was computed 
here (2, 6) and (2, 10). The absorbed insolation in the 


layer (6,10) was then computed by 
Poorman) = AG, 10) = cA(Z, 6). (4-115) 


meets noted here that the water-vapor mass above level 2 was 
assumed negligible in this study. 

Secubtiaetancen( 2.) VOj from: (A) Seeserdirect transmis 
sion of F(A) insolation impinging at the earth's surface was 
determined. The transmission of F(A) insolation is then fur- 
ther reduced by the transmissivity (1-a(G)), after surface- 


reflectance which leads to the earth-absorbed result 


<a 05 


CoOCOe Oe — ECA 2700 (2,10)Sec =z) r Caste! (1G) 


COMBE (ce) 


ivemeransmitted energy impinging upon the earth just prior to 
absorption is 

eva) sla lO (050)  [t2(G),).. (4-16 (b)) 

Meee Cloudy Cases 

In order to compute meaningful dispositions of F(A) 
meee rhation in cloudy-sky cases, cloud reflectivities and 
Sroua absorptivities after C. D. Rodgers [1967] were utilized. 
The reflectance-values used here are different from those 
Suggested by Rodgers for the F(S) wavelengths. The cloud 
memlectivities used here are RA(1) = .46 and RA(2) = .50. 
in this case there are also cloud-absorptivities to be con- 
Bewered. These were taken as A(1) = .20 and A(2) = .30, 
meopectively, In the following discussions the cloud con- 


ditions are considered totally overcast in each of layers 
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Madacated when the notations ((1, 1), (0, I), and (1, 0)) 
are utilized. | 

Semciluulesrepresemeations Of the GOmputations per- 
mommed in the various cases (({1, 1), (0, 1), (1, 0)) are 
displayed in Figures 7, 8, and 9. These figures suggest the 
parameters required in computing the disposition of incoming 
F(A)-insolation from level k = 2 to the earth's surface 
w= 10). 

Laine s/ “aopmetsSe the TdmaspOSition Ofm@a( je insoelat von 
in the overcast case with both high and low cloud layers 
mit). Figure 8 shows the case (1, 0) with an upper over- 
meeteonly. Figure 9 displays the case (0, 1) for a low 
overcast. The equations (after Warner, 1974) which relate 
ommeme parameters in these figures are listed in Appendix 
to~ee-2, B-35, respectively. 

aebiescnw ype oomeand 9, “the finsolations, AZ4, Ado, 
Peeeme.o9, and A910, etc., represent the insolation absorbed 
mimeeine layers concerned. Symbols F2, F4, F6, F8, and F9, 
etc:, depict the streams of insolation passing through the 
miameateda level, A vertical arrow implies the direction of 
insolation passage, i.e., + denotes downward insolation, + 
Moevara-reflected insolation, and ++ downward-reflected inso- 
Merron. ithe absorption quantity A{6, 8)++, for example, in- 
dicates absorbed energy remaining in (6, 8) from a downward 
ieatlected beam. 

Dicom hem EomiCctlonsuOGccun between the Garth’ s 
surface and a cloud base, or between cloud-layers in this 


model, the insolational amount reaching the lowermost reflecting 
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Figure 7 Schematic representation of F(A) insolation 


disposition in the case of two overcast layers 
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Figure & schematic representation at F(A) insolation dis- 


position with an upper overcast layer only. 
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Figure 9 Schematic representation of F(A) insolation 


disposition with a lower overcast layer only. 
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surface is affected by the downward reflected beams. Addi- 
tional insolational amounts remaining after more than two 
reflections are negligible and are not included in the models 
(Appendix B). Also note that insolation reflected upward 
meomea Lower interface, either cloud-top or ground, to the 
base of a higher cloud deck was not subjected to absorption 
BeeOr transmission through the cloud. This simplification 
resulted in slightly reduced F(A) insolation-reflectance to 
Bipace. 

Beni se rOotnemeOommote 1S that there are sthrec econ 
muewerons to the absorption of solar insolation within a 
Wayer between any two reflecting surfaces. In the Figs. 7, 
[eeeecnese Contributions are identified by arrows which in- 
miearte the portion of path being crossed. 

For each of the overcast-cloud combinations, the 
Eransmitted F(A) insolation arriving at the earth's surface 


may be defined using the following notation (TRANA) as 
TRANA(1, 1) = F1l0¥ + F10¥4¥ (4-17) 
TRANA(1, 0) = F1lO¥ + F10¥4 (4-18) 


TRANA(O, 1) 


F104 + F1l0V¥ (re Is) ) 
where the right side parameters of (4-17), (4-18), and (4-19) 
are derived in Appendices (B-1), (B-2), and (B-3). 

The F(A) insolation absorbed by the earth in each 


mase 1s derived by 


Cle) = ELON. (toa(G)) + E1014 (4-20) 
IA10(1, 0) = F104 (1-a(G)) + F10+¥ (4-21) 
WAOUneeeS 101 9 (B=o(G) le PLOY Aen)? ) 
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iaeeche last three equations, the quantity F10++ is small 
mmoueh in cach case, that no further reflections from the 


earth were considered. 


DeCOMpOsiucle (Ale Layer-Absorptions and Surface- 
Absorption Insolation 


As has been previously discussed, the standard grid- 
area weighting scheme of this study was applied to obtain 
composite values of the absorbed F(A) insolation in key 
Mayers and also within the earth's surface. The weighting 
Jaetors applied to the corresponding overcast-combination 


absorption quantities provided the following composite 


Zesults: 
Ao s(t 2 ) i= PES ee Oi ss a gp WU: 0) 
+ AZO (9 ioe eee 1s) (4-23) 
=f A264. 1) WO, ip) 
AOTC. CL(2)) = ROTO 9 gp WOO, 0) 
+ AG1O (4 9) "C1, 0) (4-24) 
+ A610 (4° Wa, Lets A610 (9 py WCe, 1) 


I 


IA10(CL(1), CL(2)) = IA10(0, 0)W(0, 0) 


SEA eIKON( te akon ( lan) 
(4-25) 
te eARONCO so.) WC Ogos 1 ) 
errame. 1yWCl, 1). 
miemwelghtang factors W(0, 0),...,W(1, 1) were first listed 
Mieequations 3-1, 3-2, 3-3, 3-4. 
4. <Absorptivity (ABA) by Layers 


Heme Uilemiitractlonal) absorptivity as well as the 


actual insolation values absorbed in the layers are considered. 


Te 





muthne Computation of absorptivity, the total undepleted in- 
Bolation at the top (k = 0) is used. The following equation 
Peo utilized in this calculation: 


FADJ = 2.00(r/r_) Cos z. (4-26) 


ABA, or the absorptivity of the troposphere was computed from 
Geeetatio of the insolation absorbed in the troposphere to 
Mmeemiisolation ineéident at the top of the atmosphere rather 


mma at k = 2: 


_ A(2, 6) + A(6, 10) ; 
nc aos 


fee bENVO (ALB) OF THE EARTH-TROPOSPHERE SYSTEM 

In considering the planetary albedo, the absorptions of 
the earth-troposphere system in both the F(A) and F(S) in- 
fomecioOnal regions must be recalled by the program. Thus 
wemderine the reflected insolational energy in F(A) as REFA, 


and it iS computed by 
REFA = F(A) - A26 - A610 - IAI0. (45210) 
Likewise REPS, the reflected part of F(S) 1s defined by 
Re = Wo.) vem > 1.0 (4-29) 


mess understood that REFA, REFS are computed for appropriate 
CL,; CL. at each gridpoint sounding considered. 


summing REFA and REFS gives REF: 
REDe = Reto REA (4-30) 


Finally the planetary albedo is related to FADJ through 


ALB = Rant ° (4-31) 
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COMPOSI Te eABS@RPTIVITY (ABG) BY THE EARTH-SURFACE; 
COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 


1. Absorptivity (ABG) of Earth 
By summing the weighted values of F(S) and F(A) por- 
tions of the incoming insolation entering the earth, the 
total insolation absorbed at the earth's surface was com- 
pereed., this quantity when divided by the extraterrestrial 
insolation gave the fractional absorptivity (ABG) of the 


earth's surface. The equation for ABG was 


LATOP +  DSHa0 


ABG = FADJ 


Ges 


where IA10, IS10, and FADJ were defined previously by Equa- 
meemss 4-25, 4-12, and 4-26, respectively. 
2. Transmissivity (ATRAN) of the Troposphere 


Also computed was the total insolational energy 


Pewee ilcident at the earth's surface just before absorption 


Ipeeene SUrface. This calculation is given by 
TRAN = TRANA + [1S10 / (1-a(G)]}. (4555) 


Here.g(G) 1s as previously defined in Equation 4-7. TRANA 

is as defined by the weighted value of TRANA(0, 0), TRANA(1, 1), 
TRANA(1, 0) and TRANA(0, 1) given by Equations 4-16, 4-17, 

4-18, and 4-19. Aiso note that the four terms of ISi0 of 

meso), (4-9), (4-10), and (4-11) have the common factor 

(1-a(G)) in the numerator and therefore each F(S) insolation 
component at the earth just before absorption need only be 
divided by (1-a(G)). TRAN may thus be viewed as the total 


msoOlational energy incident at a pyrheliometer located at 


earth. The (fractional) transmissivity of the troposphere 
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(ATRAN) is then computed from 
ATRAN = TRAN / FADJ. (4-34) 


Meee tinally that the major dispositions of the total in- 
Solation at the indicated map times have been identified by 
the fractional values, ALB, ABA, ABG, and ATRAN, each of 
Women 1S a fractional value representing the reflectivity 
falibedo), absorptivity, or atmospheric transmissivity as 
the case may be. 
3. Computational Check 

The computational scheme utilized in this model was 
checked by summing the fractional values ALB, ABA, and ABG 
mmative to the troposphere at each gridpoint. The value in 
Elemecase must total .96, since as previously noted the at- 
tenuation of solar insolation was taken as four percent as 


iemotrreamed through the stratosphere. 


Meee olATISTICAL ANALYSIS 

In order to substantiate some of the computations ere 
formed in this section, several of the most important items 
computed were statistically analyzed using linear regression 
computer programs from the BMD set of statistical programs 
mmxon, 1973). 

ie 6 Cc lear Sky Cases 

Wenger COs Ul 2). which 1s the Clear sky case of 


albedo, as the predictand, and log,,sec z and its square as 


mie predictor, the following best-fit equation resulted 


ALB(O, 0, 2) = .05705 + .54529 log,,Sec z 
(4-35) 
= .1888S5(log Sec z)?*. 


(ee 





Reese a3 
m 


where RO represents the multiple correlation coefficient. 
This result was as expected over the ocean where 
a(G) and a(R) involved logarithmic dependence on Sec z. 


Average values of ALB and log, ,sec Zz were 
NRT S20 Se log,,sec z= .34049. 


Mimer Clear Sky regression tests made use of the water-vapor 


mass path length (M) defined by 
M = (U Sec z)* log,,(U Sec z)?. (Me G8s 


This parameterization of water vapor mass is Similar to that 
developed by Hanson [1971], who used a similar M as a pre- 
Gteror in his empirical formulations of ABA for both clear 
and partly cloudy sky cases. In this study the two regres- 
mimo attempted using M as a predictor were with ABA(0, 0, M) 
and ATRAN(O, 0, M) as predictands. The best-fit equations 


mesulting were 


ABA(O, 0, M) = .11478 + .04012 M (4-37) 
R= .9675 
m 
ATRAN(0, 0, M) = .77794 - .06337 M (4-38) 
R= «8575. 


The means of ABA(0O, 0, M) and ATRAN(0O, 0, M) for the original 


93 soundings were 


ABA(O, 0, M) 


~15154 


ATRAN(0, 0, M .71987 


M 2 


i 


1.83259 (gm cm )2. 
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Jace uomsnips Between Albedo, Atmospheric Absorptivity, 


MeMooone ale I riansmissSivity. and Ground Insolation in 


the Cloudy and Clear Sky Cases 


In addition to the clear sky value of ALB, it was 
also possible to compute the composite cloudy sky case. A 
regression was formed showing the relationship between the 
mateo ALB(CL(1),CL(2))/ALB(0, 0) as the predictand and the 


meal Opaque Cloud cover (CL) as the predictor. CL in this 
model is specified by Equation 4-39 below: 


Ge (Gli(l) e+ Cl (2).= GhC1)el (2). (4-39) 


Weeeit was felt, gave a good approximation of the effective 
cloud-cover by the two layers of cloud amounts CL, and CL. 
Meeegquation 2-14 used in this model. As explained earlier 
fetnis study, a one-third reduction of both CL(1) and CL(2) 
aiter initial determination by Equation 2-14 was also tried 
in all computations involving cloud amounts. This test 
feweeesented an initial attempt at tuning the cloud model for 
mratative calculations as will later be substantiated, and 
gave results MimenOcoumatrcemememwienmtne latest litematine 
@omsatvellite reflectances. Thus all subsequent regressions 
meolving CL will also include those best-fit equations for- 
fematced using a Cl, based also upon a one-third reduction of 
Mime) and CL(Z2). A small letter subscript ‘a' in the equa- 
pron number will indicate those computations using full-cloud 
amounts while a small letter subscript 'b' will indicate those 
Metlizing a one-third reduction in the cloud model amounts. 
The first regression equations tested were those for 


ALB(CL) and ABA(CL) for the respective cases of "full" and 
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mewo-thirds cloud cover. The results are given in (4-40), 
4-41). The symbol ce ONEe again sSipenifies the multiple 


Correlation of the statistical regressions. 


ALB, (CL(1)CL(2) = ALB(0,0){1 + 3.59683CL 
= 1. 7006CL? | (4-40 (a)) 
R= .8755 
m 


ALB, (CL(1),CL(2) = ALB(O, 0)[1 + 4.13667CL 


= Seo ZC] (4-40 (b)) 
R, = +8876. 


The sample means of the values in Equations 4-40a and 4-40b 


are 
Equation 4-40a Equation 4-40b 
pipp(cei(1),CL(2)) = .4298 AUB ings Chi) = 25646 
PaeB(O, 0) = ,2035 faleB (00) Ss Al) SS 
On = 6951 Ch = S176 


The correlation coefficient in Equation 4-40 appears to show 
Strong dependence of ALB upon CL. Also note that ALB(CL(1}, 
Mee) ) values computed here are slightly higher than reported 
Dy Raschke et at., [1973], (ALB = .30) for essentially the 
femieeset Of pridpoints during the NIMBUS III period 21 January 
enrough 3 February 1970. However, the one-third reduction-CL 
Mmepears to give much closer albedo results, suggesting that 
meme full cloud amounts of Smagorinsky by 2-14 are too high. 
PECCeondmne Teco POMmconputed was between the ratio of 
meopospheric absorptivity of the cloudy air to that for the 
muucar Sky Case, versus CL. The following best-fit equation 


resulted: 
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ABA, (CL(1) ,CL(2)) = ABA(O, 0)[1 + .5043CL 


7 (4-41 (a)) 
alia? Gielea| 
R_ = .8989 
m 
ABA, (CL(1) ,CL(2)) = ABA(O, 0)[{1 + .5522CL 
Co) 
> A SCG | 
R., Ste (sy) 
Mean Statistics for this regression case are 
Equation 4-41(a) EG tase ihoimee 4. 0b.) 
ABA, (CL(1) ,CL(2)) = .1916 ABA, Pty 3c (2) = eae 
ABA(0, 0) =o .honlS ABA(0, 0) = ,1515 
CL, = .6951 CL, = soles 


In (4-41), it 1s seen that the model specifies an increase in 
atmospheric isilary) wabsSorptrvVity with increasing the cloud 
Sover, iiIhis result 1s in agreement with studies of Plante 
[1972], Warner [1974], and with the earlier study of London 
957). 

An analogous statistical regression was then developed 
momethe cloudy-sky transmissivity relative to the clear-sky 


Mm@msMissivity. The resulting regressions were 


iH 


ATRAN, (CL(1),CL(2)) = ATRAN(0,0) [1-.6020CL 


(4-42 (a)) 
LSC) 
R. = .9895 
m 
ATRAN, (CL(1),CL(2)) = ATRAN(0,0) [1-.6833CL 
(4-42 (b)) 
+ ,1246CL?7] 
Rn = .9862. 


fe 





“ 


Equation 4-42a Equation 4-42b 


ATRAN CL(1),CL(2) = .4440 ATRANCL(I),CL(2) = .4957 
ATRAN(0, 0) = .7199 ATRAN(O, 0) = ,7199 


meeEquation 4-42b is terminated at the first power of CL, the 


mesult reduces to 


ATRAN(CL(1),CL(2)) = ATRAN(0, 0)[1-.S961CL] (4-43) 


RO ao Oe 


The result (4-43) is in reasonable agreement with the Savino- 


Seestrom result which is 
Pere eeint 1), Cl(Z)) = ATRAN(O, 0) [t-€l-k) CL] ~ (4-44) 


mere kK 1S a Slowly increasing function of latitude, ranging 
meme. 55 to 0.50, according to Budyko [1956]. 
eeeeeoome Conciusions 

ihemwsctattstical analyses pertormed seem to show agree- 
ment with observational results of other investigators. Thus 
the radiation model, with the solar seasonal effects imposed 
appears to be realistic. Further, the one-third reduction in 
amounts given by the cloud parameterization model appear to 


give even better agreement with the latest satellite findings. 


pee BEDO COMPARISONS WITH PUBLISHED RESULTS 

The tropospheric albedo computations of the solar-inso- 
Mati1on model of this section have been presented for both 
ime tubl-cloud and two-thirds CL cases, respectively. These 
m@ipedos have been interpolated to whole multiples of 5° 
latitude between 20°S and 65°N. These computations were made 


mOr Gach of the four oceanic meridians and the resulting 
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zonally-averaged albedos are presented as a function of 
Haertude in fable IVY for both the full-CL and two-thirds CL 
cases, where comparison is also made with satellite clima- 
mology of Raschke et al. [1973], pertaining to the NIMBUS 
mel averaged albedo measurements over the period 21 January - 
5 February 1970. For sake of consistency the Raschke albedo 
analyses have been derived by interpolation to latitude using 
bEne same four meridians before applying the zonal-averaging 
meacess, 

The model results under the headings "FULL-CL" and "2/3- 
Om correspond to the composite cloud amounts previously 


Merced in [Table II(a) according to the equation 


= : * = 
CL CL, 43 CL, CL, CL, eso ily 


with either (CL, , CL.) Dyebauattwen 2-14 ors (Cl. ; CL5) by 
pomation 2-15. | 
iieawooesunustal results fo be observed in Table IV are 
the smali albedo-values (~ .2) between 20°S to 20°N after 
mroenke, aS contrasted with the corresponding values (in the 
mamee .5 to .4) by 2/35-CL cloud model , and more markedly with 
mmese OL the FULL-CL model (values in the range .35 to .5). 


The “hemispheric"™ cosine weighted mean albedos are as follows: 


FULL-CL, ALB = .46 
2/3-CL, ALB = .38 
RASCHKE-ALB = .30 


memeevery latitude (¢ < 50°) under comparison the 2/3-CL 
albedo computation is closer to Raschke's observations than 


mene tull-CL value. At latitudes ¢ > 55°, there is evidence 


81 





ALBEDO ALBEDO ALBEDO  FULL-CL (2/3-CL) 

Lat (FULL-CL) (2/3-CL) (Raschke) AMOUNT AMOUNT 
208 5462 .4458 | .2195 932 758 
15 5380 .4336 .2050 927 733 
10 4901 3859 .1875 867 634 

5 4334 3432  .1967 746 536 
0 4447 3478 .2325 790 556 

5 4290 3362 .2225 766 527 
10 3532 2859 . 2050 557 376 
15 3517 .2895  .2075 532 342 
20 3390 2868  .2450 443 279 
25 3664 3151 =. 2925 463 320 
30 4217 3668 .2725 559 ol 
35 4145 3703 3375 516 382 
40. 4595 4275  .3925 628 529 
45 5558 5151,  —-. 4250 833 740 
50 5682 5186 .5325 pla 759 
55 Binge "S248 5929 931 749 
60 5180 .4848 5600 764 528 
65N 4294 .4153 —-.. 7000 6 0 
pe 4622 .3820 .3010 688 514 
Avg ; 


Table IV. Comparison of planetary albedo as found by this 
memeay fOr both full and two-thirds cloud models for 16 
January 1974 and by Raschke et al. (1973) based upon NIMBUS 
III measurements. Also included are the globally weighted 
mean values at the base of the columns and the composite 


Soud amounts for both cloud cases. 
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Geeadded surface-ice reflections contributing to the Raschke 
Mealues of planetary albedo. This is especially true at the 
latitude 65°N where Raschke's results indicate ALB(65°N) = 
my, whereas the model-computations with essentially zero 
Cloud cover is close to .42. Part of the discrepancy with 
the model at high latitudes lies in the non-inclusion of a 
pemaace-1ce reflectance parameter. 

loesummarize, apart from the region covered by Sea-ice, 
there is a marked tendency for the 2/3-CL cloud model to give 
reduced albedos as suggested by Raschke et al. [1973], Von der 
Semana Hanson [1969], and others. This result of the 2/3-CL 
model is especially more clearly defined in tropical latitudes 
where Von der Haar and Hanson [1969] have discussed the re- 
eeerty Of both reduced cloud-covers and the resultant reduced 
global-albedo. These latter effects than must give rise to 
Mareecr net radiation into the oceanic surface, which they also 
Meeseribe. 

imimocctioneVi, heat-budget computations along each of 
the four meridians are presented. These computations make 
mse only ae Giheme, oClLmc OldEniOdeteLor pil poses Of radiative 
@ameulations. This selection (2/3-CL) is supported by the 
[mets Of Section VII. However, it is felt that the grid- 
pount me gaducied here for 2/3-CL is not the final answer in 
Meme roplcs. In the tropics the peculiar reflection-charac- 
MemrsticsS Of Cumulus convective elements must be more real- 
istically accounted for than by the large-scale layering 


Sonsidered in the solar insolation model of this section. 
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Peon NOo IBLE AND LATENT HEAT TRANSPORT AT THE SEA-AIR INTERFACE 


A. GENERAL 

In order to compute heat balances involving the atmosphere 
and the ocean-surface, it was necessary to compute vertical 
fluxes of sensible and of latent heat across the air-sea in- 
mieace. lhe model adopted to describe the turbulent-fluxes 
in the planetary boundary layer was essentially that used in 
the FNWC primitive equation model as delineated by Kesel and 
Winninghoff [1972] and Kaitala [1974]. The model as used in 
this study was adapted in such a way that it gave similar 
quantitative bulk-transfer results across the air-sea inter- 
face as the model used by Kaitala [1974]. This was a major 
deviation from Warner's study [1974], where he required a 
pOSitive sea-to-air sensible heat transport based primarily 


moon a Semi-empirical Bowen ratio. 


feeoeNSTIBLE HEAT TRANSPORT 

meepreviously mentioned a slightly different modelling as 
compared to that of Kesel and Winninghoff [1972] was used in 
mise study for the vertical distribution of sensible-heat con- 
_ vergence Mim eiemlaenEouOneOr lLOUU Mb in the presently opera- 
tional heat-package of FNWC, as described by Kaitala [1974] 
and by Kesel and Winninghoff [1972], the planetary boundary 


meyet sensible heat transport is given by 


K* Ty ee 
Hp (900) = 4 9C, eo) te? 
Lot as s({ 9 gt 
eG (5-1) 
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In (5-1), the parameters are as described below: 


C, = .239 cal gm ‘(deg. K)”’ 
one = Tf 
x X 
Nee = Guicwertitca mecountergnadahent lapse rate) = 


3x 10 °(deg. K)cm? 
K* = (the neutral value of the eddy coefficient) = 


=| 
PaO eme Sec 


aoa scurbulent transfer coefficient) = 
5 x 10%cm(deg. k)~* 
Pig = (surface density) = Sep Re) 


ire other parameters in (5-1) are 8g» which is the potential 
memperature at 900 mb, Oy Vitis  chemnNeOcenelalntemperacune 
meeene top of the constant flux layer, and O which is the 
Peeential temperature at the surface. 


Equation 5-1 was given by Kaitala [1974] as the planetary 


boundary layer equivalent of the bulk formula for Hp > given 
mmecrms Of Surface-layer parameters as 

Hy = 2 os ap oe Moos (5-2) 
where Cy, JS) Ghiewarac COCmErcTent , We = surface wind, and 


(T - i) is proportional to the temperature lapse for the 
mimmelayer between the sea surface and the top of the con- 
Seant-flux layer. If Hy Gimme S=Z) sare always equal, then 


mer tOllows that iS may be solved from the identity, giving 
Q 


2%, 
Ce ee ee Y 
a (Sas tan) 
K 
bz * “ps 
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where 


K = ——k-_____ 7 (5-3 (b)) 


9-10 10 
9 


An 249 7 -< = (an 10) (355° (e)5 


imaierrequirement of identical values of Hy at level x and at 


fevei kK = 9, with H. decreasing to zero at k = 8, was inter- 


i 
mmetea by Kesel and Winninghoff [1972] as¥giving rise to 
femsible-heat convergence Hes Paver. com emed at k= 9. 
meeuaely Fig, l0a supgests that the vertical convergence of 
sensible heat flux of amount Hp occurs in the layer (8, 9) 
fiem the assumption of constant heat flux in the layer (x, 9) 
is made. However, in the sensible-heat model visualized in 


mmm thesis, a linear decrease in H. with pressure was as- 


If 
moecdeabove level 'x', Fig. 10b. Note that level 'x' is as- 
Sumed to occur so close to k = 10 that one may solve (5-1) 


fOr 


H.(k = 9) = 5H, (Bulk) . - cam 


memution Of (5-4) by use of (5-1), _($5-2Z) gives the result 


2K 9,4 E 
C,V_.T. + —— - 2K y¥Y 
ae = Se ae oe ; (5-5) 
Zak 
oz 
memration 5-5 has the same format as (5-3). As was mentioned, | 


Kaitala [1974] has suggested the use K* = 2 x 105cm2sec 
mr (5-3), and this would lead to identical values as the use 
of K* = 105cm*sec in (5—-5) [hevlatter K* value was tested 


by Warner [1974] who felt, with Kaitala, that either k* value 
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was only an approximate choice of the eddy coefficient re- 
quiring at least some tuning, with a theoretical dependence 
upon the larger-scale synoptic parameters over the ocean. 

fie spavycrcal implieations of condition (5-4) along with 
the linear decrease of H,(P) with pressure to zero at level 
he=- 8 1S depicted in Fig. 10b. This condition is compatible 
Mieieuntt£orm convergence of the sensible-heat flux Hp Se the 
ameace layer into the overlying layer (8, 10). Having uni- 
mnaieconversence from level x (near surface) to level 8 is 
equivalent to maintaining a "well mixed" boundary aver Or 
mrmcKness from the surface to 8. 

With either case, Equation 5-3 or 5-5, the sensible heat 
convergences of Fig. 10-a and 10-b were considered applicable 
peeve! k = 9 and were considered to be identical. In this 
thesis the value of 2 K in (5-5) was based upon an arbitrary 
choice K* = 10°cm2sec.’ in (5-3b), rather than 2 x 10°5cm?sec ° 


In order to evaluate H Pee vemloce sot GOnGcrinyc. —ne 


T? 
surface winds Ne by use of (2-21) at all gridpoints of the 
four oceanic meridians under study. This was done using 
Pamwations 2-16 and 2-20 with contour gradients extracted from 
Peed mb field except in the latitude band +5°. Here, the 
January climatological values of gradient winds were taken 
fimom che [Tropical Wind Atlas of Atkinson [1970]. 

Besides having computed the applicable surface winds Ve 


/[meruse inl (5-2), a modification to C,. = 1.4 x 1” * (ade te tr 


D 
Meatler and Burling, 1967) was introduced, as compared with 


Cy = 2. X ee mere Sugpested by Kaitala [1974]. 
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Acter determination of wn ve @tlae holes 57, Hp was then 
momputed at each gridpoint using Equation 5-2. The values 
of Hy auceusea Pater in calculations Gt heat budgets ian 


Beerions VI and VIII. 


fe EVAPORATIVE HEAT TRANSPORT 

noma the case of Sensible heaty convergence, the latent 
heat transport by turbulence was assumed to be subject to 
unatorm heat flux convergence throughout the layer 800-1000 mb 
Meee, fib). The total amount of latent heat removed by evapo- 
ration from the ocean surface was similar to that inerr by 
Kaitala [1974] using bulk transfer theory. While a turbulent 
mem@el Similar to Fig. 1lb can be adopted, it is necessary to 
recognize that the realized latent heat may occur at arbitrary 
Mevels (for example, higher than k = 8) and therefore is not 
@ntroduced at level k = 9. The basis of the turbulent latent- 


heat flux model adopted here is shown in Fig. lib, where 


Eggo) = 22 (Bulk) (5-6) 


ood where the surface-layer DU@letranster Ormeabent heat is 
maven by 

iS E(Bulk)=Lp,_ Cy V4, 979,) - (5-7) 
From the planetary boundary-layer transfer theory analog to 
H, (900), E( 900) may be written 
= (4,.-dg) 


a (5-8) 
2g - 24 | 


Zoom) 9 Sig 
mmeto-7), and (5-8), Lis the latent heat of vaporization 


-] “ Ps 
ea aeo One Oe -601T, > Seigl fue Ty T49 aeepe Cre 
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Also in (5-8) Ag and q, are mixing ratios at the level k = 9 
mime: xX, tne top of the constant flux layer. The factor 4 
which appears in the right side of (5-6) gives the type of 
latent heat convergence as shown in Fig. 1l1b, as against 
meimad Of Kaitala (Fig. lla). Combining Equations 5-6, 5-7, 
and 5-8 leads to the solution for qd,» namely 
Cy Vs 449 + 2 az Ag 

C 


q, = (5-9) 


K 
D Ve + 2 oe 


feeeme OZ piven by (5-3(c)). As in the case of t in Equation 
mes), the factor of 2 K of (5-9) was taken equal to Kaitala's 
K of Pomeattom 5-5 (hb). Enis was done in this study by using 

a value of K* = 1. x 105cm2sec ' as against 2 x 10°cm’sec ” 
Suggested by Kaitala. Thus amounts computed for E of Equa- 
tion (5-7) were equivalent to amounts computed using Kaitala's 
eememe (Equation 5-8) only the turbulent convergence of E 

was distributed differently in the planetary boundary layer. 
Thus after computing qx (Equation 5-9), E was calculated at 


Stemecridapoint except those over ice by using Equation (5-7). 


Pree |URBULENT HEAT TRANSPORTS OVER AN ICE-COVERED OCEAN 

Any gridpoint whose surface temperature was less than -5°C 
was considered to be over an ice-covered ocean for the purposes 
of calculating H,, and Eo Im this study there were four such 
metapoints. In the handling of Hp, a formula adopted from 
Kaitala [1974] was first utilized to describe heat Donner sone 
mio thesueo (HIGE) ~from the warmer water beneath the ice 


Surface: 


HICE = BBB*(T, - T,) (5-10) 


shi 





mmere bEB 1S an ice Conduction coefficient (697.835 ergs 

cm “sec ' (deg. K)-), T, is the freezing point of sea water 
m7i.5 deg. K), and r, is the surface temperature. MHICE is 
then utilized in the following surface heat-balance equation 


memaetermine the summation of the sensible heat flux Hy plus 


the latent heat flux: 


| 2 _ re iz 
Hy + E FPeCE OA BG Fig : (Secloly 


A Bowen ratio r as described by Kaitala [1974] was then cal- 
@meated Using Equation 5-12, where the right side is an em- 
pirical least-squares best-fit applicable over continental 
amid ice-covered ocean areas. 
ote | | 

1c) a ce a DO Soins Oe SSE e Seo (aa dl7 } 
Here » is the latitude of the gridpoint, and r was never al- 
Mowed to exceed 1.33. This restraint on r meant that E was 


always bounded by 


BSS (5-13) 


(cf. Kaitala, 1974). From the definition of the Bowen ratio 
jejeas a function of Sin ¢, and using the following equation 


for the combined interface fluxes Hy + heels 


1 
Bie Eee tre (lgt =) (5-14) 


one obtains simply 


eee 
Hp = =ez (Hp + £E). (Gas) 





pemce both the Bowen ratio (r) and Hy, + E are known at each 
Emidpoint, Hp, is easily determined from (5-14). Then a simple 


Meathrangement o£ terms in the Bowen ratio definition permitted 


ee | 





tiewcempuctation of E, This is given by 


2p 


E = Te) ° (S-16) 


Equations 5-10 and 5-11 always gave negative or downward 
memes of sensible and evaporative heat fluxes at the four 
Mmenwapoints involving ice-covered oceans since QABG - Fig* 

was dominantly negative. Further, evaporation never exceeded 
so OL Hy in magnitude. These negative values appear to agree 
Pemerally with those arrived at in the work of Vowinckel and 
Maytor {|1964| which dealt in part with evaporative and sen- 


Sible heat fluxes over an ice-covered Polar Ocean areas in 


winter. 


Pome PARGE-SCALE TURBULENT HEAT FLUX ACROSS THE AIR-SEA INTERFACE 
At each gridpoint, it is possible through the adaptation 

of the turbulent boundary layer model to include heating 

rates across the air-sea interface. For use in Section VI, 


Memoeempuce the values of H. as the convergence of sensible- 


if 
Heating rate for the layer (8, 10) through ice hOM Sano 2 ent. 
Bee AlSo for Section VI, the combination E + Hp considered 
Pammerily as a heat loss at the ocean surface is computed at 
faemecridpoint by use of Equations 5-7, 5-9. For the over- 
ice cases Hp tees determined by (5-1) while Hp results 
meome(o-12) and (5-14). In any case, the model governing 


meese heating rates, H, and E has been extracted from the FNWC 


Il 
operational turbulent-boundary layer model and appended to the 
Gmeaiational model (with the 2/3-CL case) for operational com- 


Pee tvons al cach gridpoint of the four meridians of interest. 


28: 





The values of Hy. turned out to be in general small and 
negative.in low latitudes reflecting computed negative values 
Of tT, - TY. However, in the data-rich mid-latitude cross 


sections where both strong surface winds and relatively warm 


Becanm temperatures coexist, H, is poSitive and larger and 


r 
dominates the net atmospheric heating in layer (8, 10). Re- 
Baedless of the inconsistency in the sign of Hp» Lie Vvalites 


of E + H, seems to be reasonable heat-loss rates at the ocean 


If 
surface at all latitudes insofar as comparison with other 
Soumees of data {Budyko, 1956) are relevant. 

Iimatly in the large-scale sense Hy are se iieduded: as a 
mimcmiave) contribution to the heat budget of the tropospheric 


@oeunn even though it is not known precisely at what level the 


ltweont heat flux E will be realized. 
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oe MeERTONONAL CROSS-SECTIONAL DEPICTION OF THE HEAT BUDGET 
Sy ol eriecOMPULATIONS 
pee =GENERAL 
ities eeneral desion of this Section 1S £0 utrlaze all com- 
moeacelonal Subroutines explained in Sections III, IV and V 
mimcomputations of a single time-step of the heating model 
maeeloped for use in the FNWC prediction model. It was de- 


meecad atter the testing of F. and of ALB with the correspond- 


2 
mieevalues of Raschke et al. [1975] that only the two-thirds 
eroua model would be displayed in the meridional cross- 
sections depicting the single time-step computations of the 
heating package of the FNWC prediction system (Figs. 12, 
apm o( bh), l4(a), 14(b), 15(4), 15(b), 1léla), and 16(b)). 
The appropriate radiative calculations with two-thirds CL 
were thus performed at each gridpoint of the four meridians 
[iiwmror presentation purposes in this section. However, the 
wa iewtype Computations were made with full-cloud cover so 
that the mean radiative balances could be compared layer-by- 
ayer and at the levels k = 2 and k = 10. Both sets of com- 
putations were zonally-averaged and used to get mean radia- 
meomemeridional cross-sections displayed in Section VII 
mers, 17, 18(a), 18{b), 19(a), and 19{b)). 

In this study, the FNWC gridpoint processed analyses for 
erat, Lo January 1974 were used at the three Pacific cross- 
sections, while that for 1200GMT, 16 January 1974, was used 


men the single Atlantic meridian. In any case, the set of 
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meridians was considered as representative of 16 January 
gata when Converted to heating-rate calculations. 

Figure 12 depicts in symbolic language the key to the 
mempucational entries in Figs. 13, 14, 15, and 16. In 
actuality, this key list represents computations made at 
mene tl; J) gsridpoint sounding from data in the form of 
Table I(b). It covers those computations proceeding from 
miemtop Of the troposphere to the ocean surface, including 
@emeible and latent heat transfers (Equations 5-2, 5-7). 
Feresimplicity and ease of climatological data comparison, 
mmece15, 14, 15, and 16 were developed by interpolating 
earapoint results to whole 5° latitude increments: The in- 
MeLeolation routine utilized was the Lagrangian cubic inter- 


Pomaeiron Scheme. It has the form 


Se (1-2) 1-3) (ee (1-3) (15 
ee Cia eiad), ade (2a 


(LS 2 eae): ile) ue 2 ) Cie) 
CEnnes (sa Cana = 2s) 


where I is determined from an equation similar to (4-3(c)), 
Q{(I) is the desired interpolated value of the heating parameter 
at a non-integral grid-value I that corresponds to a whole 
mueetude of 5°. In (6-1), the indices 1, 2, 3, and 4 are 
moreany four increaSing succesSive integral values of I along 
gridpoint lines. The symbol I in (6-1) 1s usually applied 
Metween points 2 and 3, but this interpolation scheme may 

also be applied with good accuracy to an I-value which is 

mom -centered in the I-array (1, 2, 3, 4). In order Om CONS 
mlcte tne Interpolation to latitude, the array is applied 
sequentially to successive sets of integral I-values along each 
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B. EXPLANATION OF TERMS IN FIGURE 12 
iWecrOss-oCCtIOnN at Level kK’ = 2 
It should be explained here that all insolational 
Memameters enumerated upon in Section IV dealt with the speci- 
mie time of day that corresponded to the hour angle h at the 
instantaneous times under consideration. The resultant in- 


@rcent Solar insolation dealt with is then by Equation 4-1 
F(2) = s(—)* Cos z. 
m 


iteorder to avoid reference to specific map times, the applic- 
pee) SOlar hour-angles were: h = 35°, 10°, 55°, and 35°, 
mermepectively, for cross-sections 1, 2, 3, 4, depicted in 

fmmes, 13, 14, 15, and 16. 

Regarding the first symbol at level kK = 2 in Fig. 492, 
mete OAVE and not F(2) is displayed. This is because averages 
over the 24-hour day (16 January 1974) are considered more 
eaminetul for the data day under consideration, and hence 
meempresented in this form in Figs. 13, 14, 15, and 16. QAVE 


is the 24 hour average of F(2) found by using a similar form 





Peeeeauation 4-1 but replacing Cos z by Cos z as in (6-2) 











e Cos Z 7 
QAVE = F(2)C eS | Mor 2) 
where 
Coomaeoslteouied Sin 6 twos 6. €os oosin Hi (6-3) 
T 


tt 


H(hour angle at sunset) arc €os[- Tan ® Tan $] 


p 
6 


latitude 


solar declination for 16 January 1974. 
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Cos z is equal to the 24-hour average cosine of the zenith 
angle. The 24-hour time-averaging period is considered con- 
Sistent with the longwave flux divergences, which change 
Slowly with the time of day. This conversion to expected 
daily averaged heat-transfer quantities is compatible with 
the determination of a heat budget for the given date (e.g., 
16 January 1974). 

In keeping with the concept of using 24-hour averages, 


Other parameters for level k = 2 are defined as follows: 








QREF = REF(t)( ¢254 ) (6-4) 
where 
REF(t) = F(2) - A26 - A610 - (IA10 + IS10). (6-5) 


Here REF(t) is the instantaneous solar-reflected insolation 
Siete PTidpoint as given in 4-30 and QREF is its 24-hour 
average, assuming that planetary albedo remains constant for 
the 24-hour period, even though it was calculated for the in- 
dicated solar time at each meridian. 

The same principle will be used with regard to all 
meet) sOlar disposition parameters in converting from time- 
dependent values at solar time t to 24-hour averaged values. 
Meisein the following Equations 6-7, 6-9 below Q26(t) be- 
comes Q26, and Q610(t) becomes Q610. Superior bars ( ) will 
not be explicitly shown in the values presented in the cross- 
peetions key, Fig. 12, but will be understood in what follows. 
Finally, the 24-hour average "'system'" balance BALT is com- 
puted from 


BALT = QAVE - (QREF + F,*) (6-6) 


he 





HOjmtnc tropopause level k = 2 at the indicated latitude. in 
eomouting this balance the value of F,* was introduced from 
Poabiacwomes- tlie Net terrestrial fluxes) Fs" such as F,* that 
moapcanetmerig, I2egwere considered to be constant. throughout 
the 24-hour period which would be approximately true if the 
mloua-GOVECrS Were representative of the period. 
eee losssocccion, in Layer (25.0) 

Rerterrame tO bie. I2 themroltowime dctinitrens apply, 

In this layer all the heat-transfers shown are assumed to 


Pamencaduauivye omly in our version of the FNWC heating model. 


Pitewaveraged radiative cooling (heating) rate is given by 


BAL26 = Q26 - F26 (6-7) 


where 


Q26 daily average absorption in layer (2,6) and is 
defined relative to A26(t) by a transformation 


similar to Equation 6-2 


F26 


POR QestelarCcOOImMe 1d temo rVven by Loudt Lon 
Si 1h oe) 


Pee Coss -section in Laver (6, 10) 


Niet mGonsSidenrinesradiative transters only in the layer 


Meee the 24-hour average radiative cooling is given by 


BAL610 = Q610 - F610 (6-8) 


Where Q610 is the sum of Q68 and Q810 in Fig. 12. BALO8 has 


been taken as one-half of BAL610. However 
BAL810 = Q810 - F810 + H. = BAL68 + Hh. (6-9) 
Mere H. is the sensible heat flux transported to (8, 10) by 


i 
l@u@ation o-2. Likewise, F610 1s the sum of F68 and F810, 


fmomnesurial tilux loss values for the respective layers. Note 
here that the parameter BAL610 remained negative at all gridpoints 


100 





meeraigures 13 through 16 when considering radiative transfers 
only. 
Mepeweross-section at Air-Séa Interface (k = EO) 
ihesheat "walance at the carth’s surmacemp5ALe) con- 


sisted of the following equation 


Bee OG. - Hay* — (Het ENE (6-10) 


Here QABG is the 24-hour average insolation absorbed by the 


surface as 





WEG =FONBG(t) (Cos z/.Cos Zz} (6-11) 


(i+ Hp) is thessum Of "ie muUunpelenm neat. Lransters 
which were computed from the bulk formulas as explained in 
Peemron V¥,B. and V.C. 

Meee LOTONAL CROSS-SECTIONS OF THE VERTICAL HEAT BUDGET 

Hen tG JANUARY 1974 

meeecs 15,614, 15,0amd 16 as explained yrepresent the 
Pewee -time step of heating computations for each of the four 
meridians Sic (iteeltSmioniGy. NOG thgtet nes [our Figures 
(13-16) are divided into a and b parts with section (a) 
Beevwane “tropical' results and (b) mid-to-high-latitude eis 
Tiberi le the results depicted in these cross-sections 
meee xiibicted as representing daily-averaged values, they are 
actually based upon heat-budget parameters at the specific 
map time QQOOGMT and/or 1200GMT on 16 January 1974. Hence, 
if the results are to be used as part of the "heat-package" 
Suiproutine of FNWC, all solar-radiative absorption and-reflec- 


tance terms must be recovered as functions of GMT e.g., 
Pere) —sOnvh * icos 2 / (eos z) | 
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Refer to Fig. 


tropileal Section, 


125°W Longitudinal cross-section, 
Values computed from data for 16 January 1974 for two-thirds case only. 
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Refer to 


125°W Longitudinal cross-section, higher latitude section. 
Values computed from data for 16 January 1974. 


Figure 13(b). 
242 Or Key. 
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Refer to Fig. 


tropical section. 


170°W Longitudinal cross-section, 
Values computed from data for 16 January 1974. 
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Refer to Fig. 


higher latitude section. 


170°W Longitude cross-section, 
Values computed from data for 16 January 1974. 


Figure 14(b). 
for key. 
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Refer to Fig. 


tropreal Sseceion. 


145°E Longitudinal cross-section, 
Values computed from data for 16 January 1974. 


Figure 15(a). 


for key. 
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Refer to 


higher latitude section. 


145°E Longitudinal cross-section, 


Figure 15(b). 


Values computed from data for 16 January 1974. 
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Values computed from data for 16 January 1974. 
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and 
OMe te O bbl Cos zu (Co suze 


mee. these solar disposition terms may then be Die asl a rarerel ati 
aemection with a one-hour time-step application of the 
micrmodynamic equation of the set of primitive equations used 
meme FNWC prediction-process. The purpose for leaving all 
moors Of Figures 13,..., 16 in the form of daily-averages 
was primarily so that these figures could be compared with 
eermatology where available of the zonally-averaged heat 
budget of the ocean-atmosphere system for January. 

imididy while only the (2/5)-CL cloud model results have 
been depicted in Figs. 13, 14, 15,16. An analogous set of 
fee puedeet cross-sections were developed for the full-CL 
nee@el. These are not presented here because of the sheer 
volume of the results. However, certain radiative terms in 
both the (2/3)-CL and full-CL cross-section sets have been 
extracted from both sets for purposes of presenting zonal 
comparisons of the pmirely radiative-heat-budge@ of the FNWE 
data as it was applicable to 16 January 1974. The compara- 
mercezonally-averaged radiative cross-sections are presented 


memos, 18 and 19 of Section VII. 
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VII. ZONALLY-AVERAGED RADIATIONAL BALANCE OF THE OCEAN- 
ATMOSPHERE SYSTEM 
A. GENERAL 

iieorder €o depict the results of the purely Tadiative 
contributions, zonally-averaged over the ocean-atmosphere 
System for this study, Figs. 18(a, b) and 19(a, b) are pre- 
semecd, These figures present respectively the zonally- 
averaged radiative cross-sections for both the 2/3-CL and 
the full-CL cloud amounts. These CEOSS -sScerlonsmolowee me 
results after zonal-averaging over the four meridians con- 
Sig@ered in this study. The results are then displayed in 
the format of Fig. 17, which gives the key to the line items 
Oimeme cross-sections. Note that Fig. 17 simply combines 
Ememesddiative-transfers Q68, Q810, and also F68 with F810 
Giro. I2 to arrive at layer-average radiational warming 
maees in the combined layer (6, 10). Also in the results 
Omments Section, turbulent-transfer heating rates into the 
amemesphere (and corresponding cooling in the underlying ocean) 
Moameemerely been omitted from the results of Figs. 13,...,16, 
in passing to the radiative composites presented here. 

In obtaining a simple zonally-averaged set of radiative 
Meameing-~rates at each latitude $¢, all walues of each parameter 
mestmed in Fig. 17 at latitude 4 in the range 20S,...,65N (by 
S° steps) were simply averaged over the meridians providing 
Eme averaged value at ¢. It should be’ noted that at o = 65N, 
eiere was only onc contribution to the average. In the Southern 


Hemisphere region 20S,...,10S, only two values (from 125W and 
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Bde SoW) Of each parameter contributed to the means listed 
iteme in Figs. 18 and 19. At all other latitudes, there were 
in general four values of each radiative parameter entering 
into the computed radiative mean cross-section. AS a re- 
ie, On€ may conclude that near the northern and southern 
boundaries of Figs. 18 and 19, the listed values of the 
zonally-averaged quantities may not be exact, but the general 
woe trend 1s probably reliable. (Fig. 20 shows this to be 
meevatid result.) 

Figures 18, 19, and 20 are useful for the discussion of 
Miewearth-atmosphere system radiative balance, and in the 
discussion of the relative merits of the two-thirds cloud- 
model as compared with the full-CL model. The full-cloud 
model used by Warner [1974], is also used in comparative 
ieeseussion of the heat budget of the earth-atmosphere system 


@imoection VIII. 


ioeeenTH-ATMOSPHERE SYSTEM RADIATIONAL BALANCE RESULTS 
mbeure 20 depicts the zonally-averaged distribution of 
wiiecee parameters, Ro, Ro> and R for full cloud amounts and 
M@eenirds cloud amounts, respectively. These parameters 
fMememdetined (after Malkus, 1962) as (1) Ro» the mean radia- 
[memo ncroy transfer rate across the top of the ocean-tropo- 
Spheric system Gereured Horas Sle dtemle - mn Fig, ly: (2) 
Ro» Micminicam radiative Cooling rate in the troposphere, i1.e., 
the sum of BAL26 and BAL610 in Fig. 17; and (3) R, the mean 


Madiative warming (cooling) rate at the earth's surface, 


meee cameoras BALIOQ am Pig, 17. The relationship between 


Lak fi 
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Mremwre 20. Radiational balance at the tropopause (Ro), at 
the ocean surface (R), and in the tropospheric column (R,)- 
Solid lines denote computations made with the 2/3-CL (sub- 
scripted "1"'"') cloud model and dashed lines correspond to the 
full-CL model, (subscripted "2'"'). All computations are for 
16 January 1974 on ly(min) x 10°. 


118 





these three radiative parameters is as follows: 


R. = R+R,. 7-2) 


Table V ies esezonalinz-averaced ‘conparicoms or ee Ro» and R 
as a function of latitude $ for this study for 16 January 
m4, both for the 2/3-CL anti the full -@leasese. Aliso an - 
cluded is the en values (signifying satellite climatology) 
meer Raschke et al. [1973], and Re and R.3 Watruess ine 
latter two symbols denote December ehiiecolosieat values 
peer Berliand [1956] corresponding to R and R. Resvee ec a a 


“Hemispheric'"’ means of any Q-value are then by the following 


cosine-weighting scheme: 


65 
: oe et te Pez, 
Wt.Avg. 65 é 
0 gee 


The tabular ese of Roy R, and Ro as a function of 
miemeude and for the two cases, 2/3-CL, and full-CL, have 
Maem votted in Fig. 20, where some interesting results of 
the radiational model become evident. 

Beto cVident from Fig2920 that the Re values obtained 
mame the 2/3-CL cloud model exceed those based upon the full- 
CL model up to » = 40N. Poleward of this latitude the two 
R.-curve cross with R. (CL) becoming slightly larger (although 
Mesgatuve) than R.(2/3-CL). (ieesamnemccicmal results hold for 
both Ry and R, and for Ra EET AE 8S, Roo Upeco Jat rtudes 
m0etoN, where also cross-overs in he graphical depictions 
Pati iieiOCilncasesSweinese resuits are related to the 
Zonal ly-avyeraged radiative budgetary cross-sections of Figs. 


Tomiie eas 1S explained in Section VII (C). 


IWS, 





Earth-Troposphere Surface Net Atmospheric—-Net 


Net Radiation, R Net Radiation, R Radiation Gain 
(k=2) s (k=10) | R. 

Lat R a7 Re ye ~ Re R3 Ra R, Ra 
20S ~-1800 .0681 .0216 1622 5 2S ao094 2 O91 3 

15 mew eOa52 .0052 ~ 1624 wilekZS tL Oe ee re OWES: 

10 -1600 .0488 -.0039 1780 auld) eo Aes te! 

5 21430 .0475 .0046 902 sh oo7 = 2 One ec 2 

0 mere .0lLoz —.0270 . L802 ~1443 mes p=. Olt lO = OG S 

muozo -.0050 —.0473 1637 pdiie)/) = ol) oem. 1 Go. 

10 Wow0 —- Ole. 0428 ~ 1645 ~ 1454 ae Oe 7 oy 
15 OLD Omen0409 =-.0/ 36 ~ 1374 1226 pT (0) ala le 6 Pe 
Zee 25> —.0Sl3 —. 1010 gly 7 LOZ2 SLOSStee 1937 =. 205) weed 
pee — 0256 -.0938 -.1058 0807 0738 = lca. 1796 
ete 0400 —.0998 -.1041 ~0475 0429 “04603 yee l4 72. 9-2 147 Lee 
wee O50 —.1497 —, 1551 OS) ~0145 —) loyéys) ta dba yesa 
Pees /) —. loos —.17/06| —.0092 —.0072 =. 016) =.1597 —.8654 ~—-21597 
e-. bO25 —. 1884 —.1835); —.0094 =.0068 = GARVASMI, « SAEZ ays) 
ere e=.20027—-.1969 } =. 02/5 —.0200 =.041/) =71788 =. 1/69 —-1826 
ee 02 5) —. 1/95 —.16/71)| —.0409 —.0300 — ll 3 sro) alls yi 
ere e972 659 | eos) —. 0003 —.0509) =o t06) —. 1077 §—.2037 
pee .2G00 —.2265 -.2265| —-.0868 -.0868 eo (em le oh 
We. 
Avg 
. e500 = OOO —. 1139 0668 -0581 -0699) =.1668 -—-.1720 -—.1437 


Table V. Zonally-averaged components Be ae R. of the earth-atmosphere 
system as a function of latitude. The averaged values as found in this 
Sel@veusing subscript “1 for the 2/3-CL case and the subscript "2" for 
the full-CL case. Corresponding values of R. reported by Raschke et al. 
(1973) are denoted by subscript '4" and those of R, and R reported after 
Berliand (1956) for December climatology denoted by subscript ''3". 
Weighted averages over the latitude range O0-65N are shown in the last 


line of the table. (All values in ly min™!). 
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C. EXPLANATION OF THE MEAN ZONAL RADIATIVE-BUDGET DISTRIBUTIONS 
teeonoulae be obServed that in lable Vy .etie R. values ob- 
tained with the 2/3-CL model are greater than those resulting 
maonmene fu1l1-CL model throughout the tropics and up to lati- 
mae 40N (approximately). Poleward of latitude 40N, the 
values of R.(2/3-CL) become smaller than those for full-CL. 
imese Comparisons are illustrated clearly in Fig. 20 by the 
Zemal depictions of the Ro values calculated by the two 
@oud-model cases. Examination of Table V for comparative 
Values of the radiative balances at the surface shows the 
Same general effect at the ocean surface, with the cross-over 
between 2/3-CL and full-CL cases (Ry and R, curves ait) as ae Cena) 
meeurring near latitude 40N. Finally, the same general effect 
fomewident in the comparative values of R,(2/3-CL) compared 
to R, (CL). That these three comparisons should follow each 
maiems systematically is perhaps to be expected in view of 


Miomeonnecting relationship 


R, a Ss Ro: 


imelewchanges resulting from substituting (2/3-CL) for the 
full-CL cases, the largest gain in low and mid-latitudes oc- 
miiesein the comparative R,-curve fOr she ewsea 5 Gla)) CASe mame 1¢ 
emairest gain occurs for the computed R,-changes, anda t= 
amleditate gain is registered for the R(2/3-CL), the net 
mimmatave flux at the ocean surface. 

The same results for R.-changes may be seen by comparing 
aos tor the 2/3-CL case with Fig. 19 for the full-CL case. 


ime such a comparison it is seen that the systematic reduction 


eee 





ef CL to 2/3-CL increases the net useful solar insolation 
fem —- 2, (Figs. 17, 18, and 19) by reducimg QREF. However, 
Meare duction aepeeved=ein OQORER@2/ 3-6) tilcmantirtre Sinai iat 
latitudes ¢ > 45N. On the other hand with the 2/3-CL case, 
the F,* Melues Smet terrestrial Voss to spacemat k = 2) as 
increased by a small but consistent amount which at higher 
Buetieuaces ¢° > 45N, exceeds the very small net “solar insolla- 
mmo gain (k = 2). 

Weeeamene tore tlUrns Out that an theveonsiderationecor 
Figs. 18 and 19 of the atmospheric-layer radiative transac- 
mrons (both solar and terrestrial) that the net radiation 
Pegemetrating to the surface is influenced primarily by QREF, 
at best in lower latitudes and up to @ = 40N. Here -A(QREF) 
emeeceds A F,* aswemmcsUulteot tie chanoc simec loud covcrsanc 
mts accounts for greater R,-values. The same effects ac- 
Gemntc tor®R-changes at the surface, while at » > 45N the 
mess -over effect holds both at k = 2, and also at k = 10. 

aie curves of RS Versus) >) 1n Fie CU seticet mMeactne 
the atmospheric participation in the radiative absorption and 
memes 1on processes, follow@the lead of -A(QREF) in lower 
latitudes and of F,* in higher latitudes, which act as the 
mererne functions in the radiative balance of the system. 

Recent satellite results especially of Von der Haar and 
Hanson [1969] as well as of Raschke et al. [1973], also draw 
attention to the reduced values of global albedo, compared 
to previous heat budget estimates (e.g., London, 1957) in 


jmroprcal and subtropical latitudes. This also shows up in 


iebie Vein the comparative values of Raschke's distribution 


Mey: 72 





of R. with the "best" R. distribution computed here for the 
fo-Cl moadel. it should be noted from Table VY that the use 
Grethne 2/5-CL model leads to R.,-values which are in closer 
agreement with Reg up to latitude 40N, than that generated 
Mmyeene sull-CL model. At high latitudes in Table VY the 
comparisons between the results of these two cloud models 
Mmaeeinconclusive in their relationship to Raschke's results. 
However, his results in some of these areas are over ice, 
which enhances his effective QREF as compared with either 
model used here. 

Finally note that comparative values after Berliand 
meso] are also included in Table V for R and Ro at 10° lati- 
tude intervals based on December global climatology. The 
most noteworthy comparisons to be made here is the transi- 
miGte trom positive to negative surface net flux R, Wine noe — 
curs between 35-40N, in reasonable agreement with the results 
interpolated either from the Ry-> R,-distribution of Table V. 
In general, the global results of Berliand reflect a greater 
continental trend with latitude than is found in the 


January 1974 study over the oceans. 
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Wet. *@ONALLY-AVERAGED TROPOSPHERIC AND OCEANIC iG@AG@BUDGETS 
FORTW6 JANUARY 1974 
fee tiie TROPOSPHERIC HEAT BUDGET OVER ICE-FREE OCEANS 
Byeawerageng over the tropospheric Column tor Gace four 
feendians displayed in Figs. 15,...,16, latitude by elatitude, 
miemzonalty-averaged tropospheric balance results (Table V1 
ama Fig, 21). Values presented in Table VI and Fig. 21 are 
Ro» CiemerOpOSMnecrlcC Fadidtive Net COOlmMecedetinedmim sectron 
meer. B., and E + Hp as formulated in Ebauatioms 5-25 55-7. =e 
Maeetcimmaranmeter represents, the combined turbulent@meat -trans- 
ier Hy + E across the ocean-air interface (See) fOr serene. 
aero... 1.580). [he heat-budséet equation for the atmospiere 
column at latitude ¢ then becomes, with the notation of 


Malkus [1962]: 


aay i Sa = A eo eect Hp). (8-1) 


ere Sa ismeilou?r=ctew Stonace-hncating taueuou the Ccolunn, 
and Oia 1s the required flux-divergence of heat necessary to 
Geane about the "balance" in the column, for the observed 
meres Of R E +H, and S.. 
a i a 

lable Vi lists the zonally-averaged values of Ro» Eee H,. 
and the resulting values of oF Sa? which are computed from 
foo). As before the subscript notation “1l" and "2" denote 
ee 2/75-Cl and full-CL computations respectively. Figure 
mimdepicts cach function Ro et Hp and os + Sa? computed 
by the two cloud-models. 

hNsenotea before an section VII.C., Ra tends to exceed 


Ro MeouNMoOst of the tropical Southern Hemisphere and in the 
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Turbulent transport 










































































Tropospheric Balance 


Radiative heating heating rate in requirements 
rate RW troposphere Ces 

Lat ta 7 ae Cen ea) Ea? COR eae ea aes 
peo —.0942 -—.0913 402 21402 ~0461 0490 

ieee —- 1073 = Ao le Eo ee 17 1S 0644 0642 

Mee 2 1292 —.134/ liso ye O59 10397 0342 

cee - 1476 0. 1492 peleliD / ape =.0229 -—.0294 
feo. —.1720)—-. 0535 70859 20839 -.0782 -.0882 0069 

ely 46e—. LoG4 207 Sie ou 0950. ~—slOiesS 

ome os —. 1882 —. 1157 IBIS als74 ee aes) = 0 50G Our 
ieee —. LOO je—. 1962 bowie er ibs Bal =,Ui2o2.! =. O3o1 

Bee 9 e051 —. 1551 ~1430 .1430 = 0507 = (Uleoi) 20209 
femee- 1745 -—-1796 le eee oy -.0488 -.0538 

ame lay 2. ay — 2 1551 Mee Go's .0526 pOs20 RU SUL 
mwomee -.1648 —. 1696 -ceue 42541 0692 0644 
alo) —. 16300. 159 / = 297 eee oe ~1374 Bo Bs .0463 
eee. L/91 —.1768 a ee aay) 9s Oso -0404 
Weel, Oo —.1/6095—.1378 51901 221908 FOmES OLS YZ 0370 
See 21308 —.1371 130 1a —.0087 i=. 007@ 

fume —- 1065 —.1077 —.2037 Oneness 20098 O17 2 BOGS 
pee - 1397 -—.1397 = 0/21 —.0721 = 6 sy Sa alles 
Wt. 
Avg 

fe ooo —. 1/20)—. 1437 = (ULL ae OO: SOS 
Hem) 
Table VI. 


Zonally-averaged components Ro (E+H) and (Q. tS) of the 
tropospheric heat budget. The averaged values as found in this study 
Beime subscript "1" for the 2/3-CL case and the subscript "2" for the 
full-CL case. Corresponding values reported after Berliand (1956) for 
Weighted averages over 


(All 


December climatology denoted by subscript "3". 
the latitude range O0-65N are shown in the last line of the table. 


values are in ly min’). 





* wae w oe 4 Fn tied 





25 
SN 
Oo 
x< 
"C20 
£ 
" (E+H,), 29 
© 
. Gre 
anys 
O 
| 2 
(E +H) 25 
10 (25) 
\t 
WW 
(Oye*S), 
yy 
a Paz} 
Rat 
“29 j eee | 
10 0 40 50 69 


20 s) 


0 20 3 
Latitude (N) 


ieure 21. Distribution of Ro» Gaclalaee heel 1020 Cee as a 


Ie 
function of latitude for the two cloud-model cases. 


(solid curves -- 2/3-CL case; dashed curves -- full-CL case) 
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Nowtenern Hemisphere tor elose to 40N, MeCUEN VEGI Gree ee CaS 
a@aceaguite small. The distribution of E + Hy aS ae Lume taomarort 
matrtude 1S independent of the cloud-model except at high 
latitudes over a frozen sea-surface which occurs near $ > 60N. 


Baeure Zl shows the identical sets of E + H, up to latitudes 


P 
mie lhe curves of ae + Sa glige alll Gye) slang ing Iau ol IS 
Micmrestaual of the right sides of Equation 8-1 for both the 
wee-CL and full-ClL cases, respectively. Note that c= ‘ Sas 
Slightly exceeds (Qs 7 S432 OVer MOSt Omerhne latvtudey mance 
momcoN as a result of the corresponding of excess R,, Over 

Ro? for the same latitude range. Note also that both 

Q 


that the values Qe. it 2 anew Chica | aan O Sle vie OV Cameo (eo 


ee So distributions reach peak values at ¢ = 40N, and 
This latter feature follows from (8-1), since Qa a Sa 
mast tollow the trend of E + Hy. where the latitudinal varia- 
mmomeot the latter function peaks sharply, relative to more 
Seeless uniform behavior of R, (9). 

The hemispheric mean value of (ao = -.0160 2y min ~ 
Peeeem column. In this study, Ora has not been computed, but 
if the (CA en ts attributed to mean storage-cooling Ss, 


of the troposphere alone, the daily storage rate corresponds 


to a temperature-change rate given by 


oT (Ri + E + Hp) | 
( It y= Al oa ae x aa min/day (8-2) 
with AP mb = 800 mb in the troposphere. This resultant mean 


@n-sttu hemisphere cooling rate (8-2) over the tropospheric 


depth of 800 mb is then 


arey 





1 


fe 


ees alee © (day an 


Q 
aaa 


averaged over the cm* tropospheric column. This seems to be 
a reasonable midwinter global value. 


The@ivarlatron in the computatwon of E + Ha in theme 3-CL 


It 
and the full-CL cases which occurs poleward of latitude 55N 


Mmemexplained in Section VIII.C. 


Beeline ZONALLY-AVERAGED HEAT BUDGET OF THE ICE-FREE OCEAN 
POumene PReSeNnt, neglece the amor Varlatrons atimmmuted 
memuacen-tO-1ce conduction rates which occur in latitudes 


poleward of 55N (see for example He Of Equations 5— uae 


e 
Peevectane thas ice-conduction term, the zonally-averaged 
heat-budget of the ocean water-mass obeys the thermodynamic 


febationship 


i, eh 2 (CE es aiel)) | (GS 5) 


mom DOthe the 2/3-CL and full-CL cases. Note that Equation 
eevee is the analog of (8-1) for the atmosphere, with - (E+H,) 
fjmemmesenting generally the cooling rate for the water mass. 

iincwterms sor the raght side of (S-Simiave been computed 
for each meridian and averaged by latitude, with the results 
Siow in Tabke VII. Then an + oS ISmcOmpuced as tre residual 
Simmtehe sum of the two functions on the right side of (8-3), 
mideshown also in columns | and 2 of Table VII. The graphi- 
waamectributcion of the results listed in thas table ane 
Sievniein Piey = 22 both @for the 2/3-€h (subscript "1") and 
Micwrule-Clh (subscript "2') cases. 

As before the primary ocean-mass heating-function is R, 


EellomEadTattve heating rate at the ocean-air interface. For 


2G 





Surface Net Turbulent transport Surface Balance 







































































Radiation, R heating rate .loss Requirements 
(k=10) 
Lat Rj R2 Ra | ~ CE+H,) 1— (EFH) 2— (ETH) af CQ) #501 (Q, $5) 2(Q tS) 3 
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30 -0475 .0429 .0463 5d ehs: 1998 dE Ame ILS) 1/10) SSO tanh 
35 Ol 2 Ola 5 34 ~2341 =22186 Vos 
pom O02 —.0072 —.0116 oo) i VAs )F AL —, [003 esol oO 
oe -.0094 —.0065 cae SAN =. 2200 moe eo 
pum. 0275 —.0200 —.0417 01 00 eer 5 A INONLE Beat lis, 7/5 
pee —--0409 —.0300 301 hs}ouL se VNC) spite. QL 
p0me=.0868 —.0663 =.0509 “1156 ss = 2024" geo = 62455 
pen —.0868 —.0868 0721 10720 -.0147 -.0147 





Wt. 

Avg 

(N. aGOSue O55. 0699 
Hem) 


= ,0640°  — 7093 =..0509 


Table VIL. Zonally-averaged components R, ~ (E+H,,) and (Qiot>) of the 
surface heat budget. The averaged values as found in this study using 
Subscript "1" for the 2/3-CL case and the subscript "2" for the full-CL 
case. Corresponding climatology values for the month of December re- 
ported after Berliand (1956) are denoted by subscript "3''. Weighted 
averages are listed in the final line for the latitude range available 


in the N. Hemisphere. (All values in ly ub 
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Becure 22. Distribution of R, — (E+H, ) and of (Oe ae 
aS a function of latitude for the two cloud-model cases. 
(solid curves -- 2/3-CL case; dashed curves -- full-CL case) 
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the 2/35-CL case, the corresponding R, is considerably larger 


i 


in low- and mid-latitudes than R EO Garena C beeeace 


93 
Mierc 1S virtually no change in the -(—E + H } protriles over 
latitude (except poleward of 60N). This means that the in- 
wermal heating rate of the ocean-mass, which 1S given by the 
fieot the two terms on the right sidé of (8-3), 1s generally 
peeater im the 2/3-CL case than in the full-CL case. Here 
an iS tiles requiredsmean Oceanre heat-f£lux divergence at 
Matitude ¢. Figure 22 supports the conclusions of Von der 
ireeand Oort [1973 )ewhich requires increased lateral trams- 
bort-divergences = in the oceans, based upon the recently 
meeervyed Satellite-values of the global albedo. The values 
feeribed here to (cee uF So] SCcim HOmGOn MObOmA Gee tiles aise i- 
ite findings. 

The "hemispheric" mean value of (eg 3, )y Comnied 


between 0 - 65N is 


7 re . eet 
CO + S), = -.0859 Ly min 


Weenmecorresponds to negligible cooling rate in the water-mass 
column, although it is five times as great as the correspond- 


ing atmospheric effect. 


fee OCEAN-TROPOSPHERIC HEAT-BUDGETS FOR ICE-COVERED OCEANS 
iecomputatien or (E+ He) Over lee-Covered Oceans 

itemde termina LOM Ot sens Hy OVEYr 1ce-covered ocean 
Was computed (following Equation A-36 of Kaitala, 1974)) 


meen the following interface conditions used in Section V: 


Ea erage ag Gein 


H. 


ope a 
fee 7 BEBCT-T.) (5-10) 


131 





Mipbatitudes where 1cé-covered oceans occur, R is negative in 
January, and Dae 1S Usually pesa tive wbUliwod sdmo~Gde x. som 
magnitude smaller than the magnitude of R. Hence the (E + 
Hp) values computed at both ¢ = 60N and 6 = 65N from the 
170W meridian are both negative (here the Bering Sea is 
ieezen on January 16, 1974). 

SOLutlom for Be os and Hp + E at > = 60 and 65 N then gave 


the following results 


ive Vili. Surface heat-transpores over ice at 1/70W 


hat i ee Ly min - E + H,. Ly min 
2/7 5-268 FULL-CL 1) Sella SAE GD 

@ = 65N .0147 ,0147 re 2 ee eee 

@ = 60N .0041 .0070 = 5 US Se SIU SS 


Mees values of E + Hy Ore the tOrevOime sadism hay are USeamiin ele 
eaientation of zoOnally-averaged heat balance, both in the 
meooosphere, and the ocean as reproduced in Tables VI and 
Remi respectively. 

ities reason for thesditterence 1n compmted wadmess ot Gk + H,., 
Mesociated with the 2/3-CL and full-CL cases at latitudes ¢ 
meoeN 1S associated primarily with Equation 5-11 and the 
Peeecenice Of ice. At $¢ > 60N, the surface radiative quantities 


Saersty the inequality 


|R <a) Ren per em Gina 


2 1 
meenougch both are negative. Furthermore, eas ise Ihe daca Sic 
mean Order Of magnitude smaller than R. Ong Ri» RESPEC rave. 


hence. wien Cl > 0 


ts 2 





aes gS Pe Pe 


r rly 


mra@ethere must be a smaller turbulent flux from air to ice 
imethe Ltull-CL case. The results of Table VIII at latitude 


60N exemplify this feature. 


2. Budgetary Equations in Both Media with Ice-Covered 


Ocean 


Fiala yet 


ice represents tie Resmi ota Molec ans 


conduction transfer rate, which leaves the ocean and enters 
Piemicmospnere Over the Ice by a process not inciuded in the 


Pimoulent exchange rate E + H Hence at ice-covered grid- 


re 
points the heat budget equations of the troposphere and 


mater-body become: 


Alr-mass over 1ce: 


Qva * 34 : Ra oe Hp) Hi ce (8-4) 
Water bey snders 1ce ; 
4g * 32> Te CE eee (S25) 


mie eadditive terms on the right side of §-4, 8-5 did not 
modify the weighted-mean heating rates of the air-column or 
the water body significantly over the values resulting from 


Feees, Zi and 22. 
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ix. CONCLUSIONS 


The present work has been primarily involved with the 
development of the radiational two-layer cloud model adapted 
from Warner's earlier thesis [1974]. The data employed was 
oceanic data from the FNWC gridpoint analyses for 16 January 
iw+e Here, it was found that the net radiative flux avail- 
able for the earth-troposphere system was best verified 
against January satellite climatology (Raschke et aZd., 1973) 
Maem the Smagorinsky cloud-amounts were arbitrarily reduced 
iyea factor of one-third. It is recognized that there were 
eome minor differences in the comparisons which could be at- 
tributed to the one day space-time average as opposed to the 
average of several weeks for saeer ines data, and years inso- 
far as surface climatological data was concerned. 

Reduction of Smagorinsky cloud amounts by one third was 
found to provide more realistic global albedo estimates than 
Meeuilted using unreduced elond aise s. In the comparison 
With satellite data, the radiative model with a reduced CL 
gave slightly higher albedo estimates in the tropics and 
mid-latitudes and slightly lower albedo in the high-latitudes. 
Sire radiative model still lacks an effective surface albedo 
model over ice. 

With the reduced cloud amounts greater surface net heating 
rates in the low to mid latitudes were computed than when 
using the original Smagorinsky cloud amounts. This latter 
imac 15 in accord with recent heat balance studies using 


Sometime data which have found an increase in surface 
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fowls particularly in the tropics. It is felt that further 
investigation of the tropical radiative balances are neces- 
sary. Thought should also be directed towards the geometry 
of solar radiation reflection problem among cumulus subgrid 
Sells in the tropics as contrasted with the larger-scale two- 
Cloud layer model used in this study, which proved fairly 
metiable in middie and high latitudes. Moreover, specifica- 
fon OL the primarily cumuliform cloud distributions in the 
mEepi1cs 15 evidently closely related to the subgsrid scale 
mem ection model rather than to the large-scale distribution 
of relative humidity (of Smagorinsky). It is recommended that 
for~thner testing of the radiative model for the different sea- 
jews Of the year continue, particularly with the one-third 
Or Similar reduction in the Smagorinsky cloud-amounts. 
berther differences from climatological data might sbe 
meaeeeputed to the treatment of cloud reflectivity and absorp- 
tivity for solar wavelengths, as well as the parameterization 
emene sca-alr exchanges of sensible and latent heat. in the 
[mieeer case a new interpretation of the simple bulk formulas 
ends of a redefined distribution of the turbulent sensible and 
fiarent fluxes in the planetary boundary layer parameterization 
model seemed to give a new, conceptually satisfying view of 
the planetary boundary layer model. Additional testing with 
the redefined planetary boundary layer or any new interpre- 
Mieoneort such, might also prove productive in providing bet- 
feonmmcommanlsons with recent investigations such as the GATE 


enaebe@MEX €xperiments involving surface heat budget calculations. 
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APPENDIX A: EMISSIVITY E. FORMULAS AS A FUNCTION OF WATER 
VAPOR AND CO.“ ABSORBER MASSES ALONG PATH OF 
INTEGRATION. 


The general form of the calculation required for computing 


met LR flux is: 


216 


Ewo (UCT) aB Cea 


5) 

I 
td 

i 
i 


0 


where Bio = Stephen-Boltzman blackbody flux at Ts and € 


Meevhe emissivity of the atmosphere as a function of both 


WC 


water vapor and CO. absorber masses along the path of 

integration. The hatched areas in Section III figures are 
Bio a 

represented by i E,. GB which is the key integral to be 
B=0 : 


Semeuted. in order to perform this calculation, Sasamori 


[1968] proposed emissivity formulas of the following form: 


SAO 38) = rte log, fU(6,8)] + pees 
+ .07262{(1. = .62556[U(6,8) + -02861°2°) 


% [log, ON Goats) Pt ale orga : Che?) 


for the layer 6-8 relative to the level 6 as reference level 


mer temperatures T 7 210 deg k where ¢€ is very nearly 


we 
temperature independent. Also note that U is the water 
vapor absorber mass given by (2-7) measured from the 


reference level 10, 6, or 2 and C is the correspondingly 
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carbon dioxide absorber mass. Water vapor emissivity alone, 
eontribution is represented by the quantity within the brace 
of Eq. A-2. The second brace indicates the additional contri- 
Bel On made by the part of the CO. emissivity transmitted 
through the water vapor "overlap" near 15um, i.e., the part 
not already counted in this spectral region by the water 
vapor dntegration. 

Equation A-2 was used to compute atmospheric emissivities 
meee the layers (8,10), (6; 10),°@)) 10).3. = Cle ake) 
hom When the running mass values U(8, 10), ¢€(8, 10), ete. 
were ready to be implied from the radiative soundings. To 
complete the Feat integration given by (A-1) from level k = 1 
fmomene Origin (B = 0, sen 1.0), on the Yamato Radiation 
Chart (Appendix Fig. A-1), U(T,); C(T,) soundings were 
adopted as isothermal from T = Ty (k = sete, T = a (k = 005 
and then extended from the last "sounding" point (U[(0, 10), T,]) 
along the last emissivity value for T < T, = 210 deg k where 
elie emissivity of the water vapor is known to be temperature- 
dependent. This value of emissivity is denoted ¢ and is 


We 


written, after Sasamori [1968] 


Eee. 10), 7) 1) = 8.34 p,£ +353 logo U(0,10) - .44) 


x u(0,10)b7 +93455 log; UCO.10) - .705] 


$1807 (2353) Togig, U09.10) 4 3.56) ) 
+ ,07262{(1. - .62556[U(0,10) + .0286}°°°) 


: [log.9 C@@B10) + 1.064 )} (A-3) 
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Figure A-1. Depiction of Ss given by Eq. (A-3) on the 
Yamato Radiation Chart. Area above heavy dark curve repre- 


sents 2 above level k=1 from a radiative sounding. 
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Mere the final term in the braces of Ea. A-3 is almost the 
same itorm of the CO, Coreecl lone Te the total emissivity is 

in Eq. A-2. In Eq. A-3, however, this term refers to the 
full layer (0, 10), the first term is temperature dependent 
as well as dependent upon the summed absorber masses relative 
to the reference level k = 10. It can be seen clearly that 
Eq. (A-3) is usable with proper final parameters U(0, k), 
eee ik), and T, for the determination of the curvilinear 
merkeion of the downward flux contribution through level k 

for all of the reference level cases to be investigated below, 
met 2S k = 2, 6, 

Itc hnould also be noted that Soa has been integrated over 
feeevemperature T < Ty, for which the water vapor emissivity 
curves are temperature dependent (see Appendix Fig. A-1) so 
that the result of Eq. A-3 is really the average of the 


final curvilinear emissivity value multiplied by BL = St t,/. 
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Masibilx Bs FORMULAS USED IN DEPICTING THE DISPOSITION OF 
INCOMING SOLAR (F(A)) FROM LEVEL k = 2 TO THE 
PARTH S SURE ACE 
These formulas are depicted by Figs. 7, 8, and 9 in 
section IV. They are used in this study to calculate amounts 
of insolation absorbed by the atmospheric layers and the 
earth's surface. For further explanation and key to these 
ioemilas consult Section IV.C. Dashed separation lines are 
used to divide those equations involving absorptions in 
mecessive layers e.g., (2,4), (4,6), (6,8), etc. 


fWeecase (1,1) Overcast in Both High and Low Cloud Layers 
Gece also Fie. 7) 


FY 


Il 


PCA) (1-.2710(2,4) See 217° > 


Fu4 = FUY(RA(1)) 


CBSE) 
AQUY = F(A).271(U(2,4) Sec 23° 393 
Agus = Pus.271[U(2,4) See 29°39? 
F6 = F44(1-RA(1)-A(1)) 
AHU6 = FUY(A(1)) (B-1b) 
AD6 = ADUY + ADUA + ANC 
TpD68 = 1-.271(U(6,8)5/31° 293 

(B-2) 


T™D910 = 1-.271[U(9,10)5/3]° 293 


oe 


where TD is a term to denote the Manabe-Moller transmissivities 


nmonwinselation beneach a cloud layer. In this event, Sec 2 


140 





is changed to the mean slant-path sec z = 5/3 [after 


Katayama, 1966]. 


F8¥ = F64(TD68) 

F8t = F84(RA(2)) 

FO¢+ = F8t(TD68) 

F6++ = F6+4(RA(1)) 

F844 = F6++(TD68) (B-3a) 
A68+¥ = F64 - F8¥ 

A68+ = F8t - F6t 

A68++ = F644 - F844 

A68 = A68Y + ABBt + ABBY4Y 


FO = F8+(1-RA(2) - A(2)) + F8+¥+¥(1-A(2)) 


A89 = F8+¥(A(2)) + F8+¥4(A(2)) a 
FlO¥ = F9¥(TD910) 
FlOt = F10+(a(G)) 
FOt = F10+(TD910) 
F9¥+ = F9O+(RA(2)) 
FlO+¥+ = F9¥4+(TD910) | (B-3c) 


A910+ = F9OY - F1OY 

AQ1O0¢ = F1lOt - FY9t 

AQ1O+YY = FOYY - FILOVY 

A910 = AQ10¥ + AQ1Ot + AQLIOYY 


A810 = A89 + A910 





Case (1,0) Upper Overcast Only (see also Fig. 8) 


Hele = FCA). e7MiuCe, 4) see zie” 

AQUA = F(A){1-.271[U(2,4) Sec 2}°3°3}RA(1) 

303 (B=44) 
xfi27 lime. 4) See Zi } 

A2u = Aoue + A2ut 


F6+ = F(A){1-.271[U(2,4) Sec 29° 2°3}f1-RA(1)-A(1)] 

ANG = F(A){1-.271[U(2,4) Sec 29° 2°3}a(1) (B=) 
A26 = A2H + ANG 

F104 = F6+(TD610) 

F10+ = F10¥(a(G)) 

F6+ = F104(TD610) 

F644 = FO6t(RA(1)) (B-Uc) 
FlO¢¢ = F64+(TD610) 

A6104 = F6+ — F104 

A610+ = F10+ - F6t 

A6104+¥ = F6¥¥ - F10+¥4 


A610 = A6104 + A610¥ + A6104¥4¥ 
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Case (0,1) Low Overcast Only (see also Fig. 9) 


A264 = F(A).271[U(2,6) Sec 29°273 

A68+ = F(A).271{([U(2,8) Sec 21° 293-fu(2,6) Sec 21° 393} 
F84 = F(A){1-.271[U(2,8) Sec 2]° 2°93} 

F8+ = F8¥[RA(2) ] 

F6t = F8+{1-.271[U(6,8) Sec 2] 293} (B-5a) 
F2t = F8t{1-.271[U(2,8) Sec 2]°273} 

A68t = F8+ - F6t 

A26+ = F6t+ - Fot 

A26 = A264 + A26+ 


Om shorii—RA(2))- AC@)] 


A89 = F84¥[A(2)] _ 
F1O¥ = F94(TD910) 
F10t+ = F10+(a(G)) 
FO+ = F10+(TD910) 
FO¥+ = FOt(RA(2)) (B-5¢) 


F10++ = FO4+4(TD910) 

A9Q10+ = FOY - F1O¥ 

A9104# = F1Ot - FY? 

AQI1O¥Y = FOYY - FIOYVY 

A910 = AQ10+¥ + AQ1Ot + AQI1OV4Y 


A610 = A68Y + AB68t + ABI + AQIO 
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SUBROUTINE FOR COMPUTATION OF SURFACE WINDS FROM SECTION II.C. 
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SUBROUTINE FOR COMPUTING TEMPERATURE AT LEVEL 
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SUBROUTINE FOR COMPUTING LATENT HEAT FLUX IN THIS STUDY FROM EQ. 
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SUBROUTINE EOR COMPUTING SENSIBLE HEAT FLUX IN THIS STUDY FROM EQ. 
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